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ABSTRACT 

A 96-ks Chandra X-ray observation of p Ophiuchus cloud core A detected 
87 sources, of which 60 were identified with counterparts at other wavelengths. 
The X-ray detections include 12 of 14 known classical T Tauri stars (CTTS) in 
the field, 15 of 17 known weak-lined TTS (WTTS), and 4 of 15 brown dwarf 
candidates. The X-ray detections are characterized by hard, heavily absorbed 
emission. Most X-ray detections have visual extinctions in the range Ay ~ 
10 - 20 mag, but several sources with visual absorptions as high as Ay ^ 40 
- 56 mag were detected. The mean photon energy of a typical source is (E) 
~ 3 keV, and more than half of the detections are variable. Prominent X-ray 
flares were detected in the unusual close binary system Oph SI, the X-ray bright 
WTTS DoAr 21, and the brown dwarf candidate GY 31 (M5.5). Time-resolved 
spectroscopic analysis of the DoAr 21 flare clearly reveals a sequence of secondary 
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flares during the decay phase which may have reheated the plasma. We flnd that 
the X-ray luminosity distributions and spectral hardnesses of CTTS and WTTS 
are similar. We also conclude that the X-ray emission of detected brown-dwarf 
candidates is less luminous than T Tauri stars, but spcctroscopically similar. 

Simultaneous multifrequency VLA observations detected 31 radio sources at 6 
cm, of which ten were also detected by Chandra. We report new radio detections 
of the optically invisible IR source WLY 2-11 and the faint Ha emission line 
star Ehas 24 (class II). We conflrm circular polarization in Oph SI and report 
a new detection of circular polarization in DoAr 21. We flnd no evidence that 
X-ray and radio luminosities are correlated in the small sample of TTS detected 
simultaneously with Chandra and the VLA. 

We describe a new non-parametric method for estimating X-ray spectral prop- 
erties from unbinned photon event lists that is applicable to both faint and bright 
X-ray sources. The method is used to generate fx, logT and Lx light curves. In 
addition, we provide a publically-available electronic database containing multi- 
wavelength data for 345 known X-ray, IR, and radio sources in the core A region. 

Subject headings: infrared: stars - ISM: clouds - ISM: individual (p Ophiuchi 
cloud) - radio continuum: stars - stars: formation - stars: pre-main-sequence - 
X-rays: stars 
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1. Introduction 

The p Ophiuchus cloud is one of the nearest sites of on-going star formation. The optical 
and infrared study of Chini (1981) gave a mean distance d = 165 ± 20 pc, but more recent 
work based on Hipparcos and Tycho data for off-core sources (Knude & Hog 1998) suggests 
the value may be as low as 120 pc. We use d — 165 pc in this paper. Despite its proximity, 
very few stars are optically revealed due to high extinction. Even so, previous observations at 
X-ray, near-infrared, far-infrared, and radio wavelengths show a diverse collection of objects 
including protostars, T Tauri stars, and brown dwarf candidates (BDCs). Near the center 
of Cloud Core A (as designated by Loren et al. (1990)) are the magnetic B4 star Oph SI 
(Andre et al. 1988, 1991), the Class radio source LFAM 5 (Leous et al. 1991), and the 
class I protostars GSS 30 IRS-1 and GSS IRS-3, first identified by Grasdalen et al. (1973). 

The first X-ray images of p Oph obtained with the Einstein Observatory revealed a large 
population (Montmerle et al. 1983) of flaring young stellar objects (YSOs). In subsequent 
observations with ROSAT, ASCA, and Chandra, Class I protostars, Class II classical T- Tauri 
stars, and Class III weak-line T-Tauri stars ^ in p Oph have exhibited some of the largest 
stellar flares ever observed (Koyama et al. 1994; Casanova et al. 1995; Grosso et al. 1997; 
Kamata et al. 1997; Grosso et al. 2000; Tsuboi et al. 2000; Grosso 2001; Imanishi, Koyama, 
& Tsuboi 2001; Imanishi et al. 2002). 

The ubiquity of X-ray flaring and rapid variability suggests that magnetic flelds play a 
central role in the X-ray emission of most young stars. However, the details of this process 
and the geometry of the magnetic structures involved are still under intense study. Evidence 
for extended magnetic structures around YSOs comes mainly from high-resolution radio 
observations. VLA observations have detected moderately circularly polarized nonthermal 
emission indicative of ordered magnetic flelds in WTTS such as Hubble 4 (e.g. Skinner 
(1993)). At higher angular resolution, VLBI observations have shown that the radio emis- 
sion of WTTS typically arises in compact regions of diameter <25 with high brightness 
temperatures T5 > 10^'^ K (e.g., Andre et al. 1992), (Phillips et al. 1991). Of particular 
interest here is the strong nonthermal radio source Oph SI in p Oph core A whose magneto- 
spheric diameter D 12Ri, has been directly measured using VLBI techniques (Andre et al. 
1991). Oph SI consists of a magnetic B4 star (Oph SI A) and a K = 8.3 companion (Oph 



Class I, II, and III designations for YSOs are based on the value of the spectral index a = 
d \og{\Fx)/d logA, which is usually evaluated in the near to mid-IR spectral range. As discussed by Haisch 
et al. (2001), class I sources have a > 0.3, flat spectrum sources have —0.3 < a < 0.3, class II sources have 
— 1.6 < a < —0.3, and class III sources have a < —1.6. By comparison, class sources are usually not 
detected in the near-IR. 
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Sl B) at a separation of 20 mas discovered by lunar occultation (Simon et al. 1995). The 
B4 star is unusual, belonging to a small class of magnetic OB stars which includes 9^ Ori C 
(07 V), the central star of the Orion Nebula. 

Although the existence of extended magnetic structures in young stars is now firmly 
estabhshed, it is not yet clear how these magnetic structures relate to the observed X-ray 
emission (and its variability) nor is it even certain that the emission processes in YSOs are 
the same as in older late-type stars like the Sun. In the Sun, X-ray emission is ultimately 
linked to the presence of dynamo-driven magnetic fields. X-ray emitting plasma is confined 
in magnetic loops and X-ray fiarcs arc an observational consequence of rapid energy release 
during magnetic reconnection events in the corona. However, in younger stars the situation 
is more complex since magnetic fields may be of primordial origin or produced by dynamo 
mechanisms that are different from the Sun, and interaction of the stellar magnetosphere 
and the surrounding disk may be possible. 

Shu et al. (1994) and Armitage & Clarke (1996) have invoked extended magnetic fields 
that couple the photosphere to the disk in order to explain the slow rotation of classical 
T Tauri stars. If this basic model is correct, then the magnetic fields must extend to the 
Keplerian co-rotation radius, several stellar radii from the photosphere. Quasi-periodic count 
rate variations in objects such as the Class I source YLW 15 in p Oph core-F (Tsuboi et al. 
2000) have been interpreted in terms of differential rotation of the protostar and its accretion 
disk, leading to magnetic shearing and reconnection thereby producing large-amplitude, 
long-duration, quasi-periodic X-ray flares (Montmcrlc et al. 2000). Also, high-resolution 
Chandra grating spectra of the CTTS TW Hya show weak forbidden-line emission that may 
be attributed either to high electron densities (Kastner et al. 2002) or to UV irradiation 
of the X-ray emitting plasma (Gagne et al. 2002), possibly from a UV hot spot. Either 
explanation suggests that material may be funneled inward along field lines from the inner 
edge of the disk to the photosphere, producing an X-ray emitting shock at or near the stellar 
surface. 

The possible effects of circumstellar disks and/or accretion on the X-ray emission of 
young stars is a subject of current interest and debate. The presence of a disk is usually 

inferred from non-photospheric excesses at near to mid-IR wavelengths. In a Chandra study 
of the Orion Nebula Cluster (ONC), it was found that stars with and without IR excesses 
have similar X-ray luminosity functions (Feigelson et al. 2002, 2003). This suggests that IR 
disks themselves have little if any effect on X-ray emission levels. 

However, it is not necessarily the case that all stars surrounded by IR disks are actively 
accreting, and there is now some evidence that accretion itself (rather than IR disks) may 
have an effect on X-ray emission. Accreting sources are usually identified by UV excesses 
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or by optical diagnostics such as emission at Ha, O I, or Ca II. Such optical diagnostics are 
not generally available in heavily- obscured regions such as p Oph, but are accessible in other 
star-forming regions. The i? 05*^ T surveys of Taurus- Auriga (Stelzer & Neuhauser 2001), 
NGC 2264 and Chameleon 1 (Flaccomio, Micela, & Sciortino 2003) and Chandra observations 
of the Orion Nebula Cluster and IC 348 (Flaccomio et al. 2003; Preibisch & Zinnecker 2002; 
Stassun et al. 2004) find that stars with high accretion rates have, on average, 2-3 times 
lower X-ray luminosity than stars with low accretion rates. Moreover Flaccomio, Micela, & 
Sciortino (2003) show that, in the Orion Nebula Cluster and in NGC 2264, stars with infrared 
excesses have lower Lx and Lx/-Lboi than stars with no IR excess, although the effect is less 
pronounced. Flaccomio et al. (2003) argue that accretion and a disk affect the magnetic 
geometry, thereby reducing the available surface area at the photosphere for closed magnetic 
loops, thus reducing Lx and Lx/Lboi- Flaccomio, Micela, & Sciortino (2003) further suggest 
that disk dissipation and/or decreased mass accretion causes an increase in Lx/Lboi during 
the first few Myrs of PMS evolution. After 10 Myr, Lx/Lboi declines with rotation. 

We present simultaneous observations of the p Oph A cloud with the Chandra ACIS-I 
camera and the VLA that provide new information on the cloud population and its X-ray 
and radio properties. Our primary objectives were (i) to obtain a sensitive high angular 
resolution X-ray and radio census of the young stellar population in p Oph A and quantify 
the X-ray and radio properties of the population, (ii) to accurately identify counterparts to 
X-ray sources at other wavelengths, particularly in the near-infrared and radio, (iii) to deter- 
mine if significant differences exist in the X-ray properties of the different stellar subgroups 
(Class O/I, II, III and brown dwarf candidates), and (iv) to search for possible relationships 
(including correlated variability) between the X-ray and radio emission of YSOs in p Oph A. 
These simultaneous high-angular resolution X-ray and radio observations are so far unique 
and provide compelling evidence for strong magnetic activity in the embedded p Oph A 
population. 



2. A Database of Young StellEir Objects in p Oph A 

As part of this study, we have compiled an electronic database consisting of 345 objects 
in the Chandra AClS-1 FOV. It includes previously cataloged infrared. X-ray, and radio data 
and visual extinction estimates, as well as new data presented here. All figures and tables 
herein are based on information in this database. The database is publically available and 
is described in more detail in Appendix A. 
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3. Chandra Observation, Data Reduction, and Source Identification 

3.1. Chandra Observation and Data Reduction 

The Chandra X-ray Observatory observed the p Ophiuchus Cloud A (L1688; Lynds 
1962) continuously for 96.5 ks, beginning at 2335 UT on 2000 May 15 and ending at 0319 
UT on May 17. Chandra collected consecutive exposures of 3.24 s on the Advanced CCD 
Imaging Spectrometer (ACIS), using ACIS chips 10-13 and S2-S3. The ACIS-S chips were 
located far off-axis and were not used in our analysis, and will not be discussed in this paper. 
The high-resolution mirror assembly (HRMA) and AClS-1 camera are described in detail in 
the Chandra Proposers' Observatory Guide ^ and in Weisskopf ct al. (2002). The ACIS 13 
aim-point was on Oph SI at (J2000) a = 16^26"^35?2, 6 = -24°23'13". 

We applied a standard data reduction using the Chandra Interactive Analysis of Obser- 
vations software ^ CIAO v2.2.1. The ACIS data were obtained in VFAINT mode, allowing 
us to remove afterglow events. The event list was filtered to include events with standard 
ASCA grades and with photon energies in the 0.5-7.0 keV band, thereby significantly reduc- 
ing the particle background. We checked for background fiares and none were found. Source 
identification (§3.2) and fight curve variabifity analysis (§5.3) are described below. Spectra 
of selected sources, along with source-specific ARF and RMF response files, were extracted 
using the CIAO tool PSEXTRACT. 

3.2. Chandra Images and Source Identification 

Two images which contain the entire 4-CCD ACIS-I FOV were generated to search for X- 
ray sources: a full- resolution 2800x2800 pixel image (0.492" pixels) and a binned 1400x1400 
pixel image (0.984" per pixel). Figure 1 shows the full-resolution 2800x2800 pixel image. 

To register the Chandra image against infrared positions, we selected 43 2MASS sources 
with obvious counterparts on the Chandra image. The 2MASS J2000 positions and Chandra 
physical pixel positions were then used to derive a 4-coefficient ACIS-I plate scale solution 
using the Starlink program ASTROM. The ASTROM — nominal offsets in R.A. and Decl. 
were —0.84" and —3.68", respectively. These offsets were applied to the Chandra image 
header keywords to obtain the coordinates fisted in Table 1. 

We detected X-ray sources on each image using CIAO's wavelet-based tool, WAVDE- 



^ htt p : / / cxc . harvard . edu / udocs /docs/ docs.html 
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TECT. We ran WAVDETECT using scaling factors 2, 4, and 8, with a false alarm probability 

of 10~^. In some cases, regions near the edge of the field of view were detected as sources by 
WAVDETECT because the off-field background is zero. After visual inspection, these spu- 
rious detections were omitted. Two X-ray sources in Table 1 were not detected by WAVDE- 
TECT: SKS 3-12=J162622.2-242447 and BKLT J162636-241902=J162636.8-241900. They 
have been added because they contain more than 7 counts and correspond to known IR 
sources. 

WAVDETECT computes the 1-a major- and minor-axis radii of each source ellipse. 
Events within a 3-a major- and minor-axis radii ellipse were extracted for each source. In a 
few cases, the extraction ellipses were shrunk to avoid region overlap with nearby sources. 
The net counts given in Table 1 are the raw counts in the elliptical extraction region minus the 
estimated background in the same region. To determine the number of background counts 
we used the area of the extraction region and assumed a constant 0.5-7 keV background rate 
over the CCD of 0.17 counts s~^ per ACIS-I CCD (see the Chandra Proposers' Observatory 
Guide). The typical background during our 96.5 ks exposure is approximately 4 counts per 
source. Table 1 lists the 87 sources in the 17' x 17' AClS-1 FOV. Of the 87 Chandra sources, 
60 are identified with catalogued IR or radio sources. All but three of these identifications 
have X-ray — IR offsets < 1". 



4. VLA Observations and Data Reduction 

Multifrequency radio observations of p Oph core A were obtained simultaneously with 
Chandra on 2000 May 16 using the NRAO^ Very Large Array. The observations provided 
continuous rise-to-set monitoring of p Oph over ~7 consecutive hours, as summarized in 
Table 2. The VLA was in C-configuration and all 26 operational antennas were used at each 
observing frequency (no subarrays). The primary frequency was 4.86 GHz (6 cm), but briefer 
coverage was also obtained at 1.42 and 8.46 GHz to search for additional radio sources that 
might have been missed at 4.86 GHz. Since none were found, the discussion below focuses 
mainly on the 4.86 GHz results, which provide the best sensitivity and most complete spatial 
coverage. 

Nine separate 6 cm pointings were distributed on a grid centered on Oph SI with a ~ 4.5' 
grid spacing in order to provide nearly uniform radio sensitivity across the full Chandra FOV. 
Data at each grid point were obtained in two ~15 minute scans, with each scan bracketed by 



^The National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation, 
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-8- 



phase calibrator observations. In addition, we obtained a single 1.4 GHz pointing centered 
on Oph SI and four 8.4 GHz pointings centered on Oph SI and the known radio sources 
LFAM 2, 5, and 15 identified in the previous VLA survey of Leous et al. (1991, hereafter 
LFAM). 

We edited and calibrated the data using the AIPS^ software package. The data at each 
grid pointing position were combined to produce separate cleaned maps at each frequency 
using the AIPS task IMAGR with natural weighting. In those cases where the field contained 
a bright radio source, we used phase-only self-calibration to improve the dynamic range. 
Typical rms noise levels at phase center in Stokes I cleaned maps were 36 /xJy/beam (6 cm), 
40 /iJy/beam (8 GHz), and 220 /xJy/beam (1.4 GHz). The higher noise level at 1.4 GHz is 
due mainly to the presence of a bright extragalactic source (BZ6) located ~ T southeast of 
Oph SI. 



4.1. VLA Source Identification 

We identified radio sources using cleaned maps at a 5(7 detection threshold, yielding the 
31 detections listed in Table 3. We measured peak fluxes with the AIPS task IMEAN, and we 
measured total fluxes using IMFIT (Gaussian source model) and TVSTAT (pixel summation 
inside of the 2(t flux contour). Fluxes for off-axis sources were corrected for primary beam 
attenuation using PBCOR. 

Using our database, we attempted to identify known counterparts for each of the radio 
detections in Table 3. Previous detections at one or more wavelengths were found for 18 of the 
31 VLA sources (17 of 28 in the Chandra FOV). Sixteen of the radio sources are conflrmations 
of previous VLA detections obtained by LFAM and Stine et al. (1988, hereafter SFAM). We 
note that the VLA beam did not resolve Oph SI A and SI B. As a result, we assign the 
radio emission to the binary Oph SI. 

We report new radio detections of the IR source WLY 2-11 (J162556. 1-243014) and 
the emission line star Ehas 24 (=WSB 31=J162624.0-241613). WLY 2-11 is an embedded 
IR source with no visible counterpart and an estimated extinction of A\ = 7.3 (Wilking 
et al. 1989). Elias 24 is a faint {R = 18.5) star whose H« emission is barely detectable 
above the continuum (Wilking et al. 1987). Its spectral type based on near-IR spectra is M2 
(Doppmann, Jaffe, & White 2003). 

Two relatively bright 6-cm radio sources observed by LFAM were not detected in our 



^Astronomical Image Processing System (AIPS) is a software package developed by NRAO. 
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VLA observations: LFAM 20 and LFAM 22 in Table 1 of LFAM, whose respective flux 
densities in April 1988 were 824 /xJy and 599 /xJy. Our upper limits of ~150 /iJy (Sa) for 
these two objects indicate that they have weakened considerably. 

Thirteen of the VLA detections have no known counterparts and were not detected in 
previous VLA surveys of comparable sensitivity (LFAM, SFAM). One of these, J162615.5- 
243428, is relatively bright with S^cm = 4.09 mJy and is offset by ~ 6'.'7 from the HII region 
Oph 11. Given that these unidentified sources were not detected previously, some may be 
variable emitters associated with deeply embedded objects. 

5. X-ray Source Properties 

5.1. Summciry of X-ray Properties 

Chandra detected X-ray sources with luminosities spanning more than four orders of 
magnitude from log Lx = 27.67 - 31.82 ergs s~^ (§5.2). The mean photon energy of detections 
is typically (E) 3 keV, but the two hardest sources have (E) 5 keV (Table 1). Variability 
is detected in 45 of 87 sources (52%) as judged by a value of the Kolmogorov-Smirnov statistic 
KS > 1 (§5.3). For those detections with known IR counterparts (Table 4), visual absorption 
ranges from Ay — 1.8 - 56.0 mag. The faintest IR source detected by Chandra was GY 144 
(K = 13.46). 

The list of X-ray detections in Tables 1 and 4 includes 12 of 14 known GTTS and 15 of 
17 known WTTS within the AGIS-I FOV. The X-ray luminosity distributions of GTTS and 
WTTS are similar (§7). Four of 15 BD candidates were detected. Apart from their lower 
X-ray luminosities, the X-ray properties of BD candidates are otherwise similar to TTS (§9). 
The class radio source LFAM 5 and the class 1 protostars GSS 30 lRS-l/lRS-3 were not 
detected by Chandra, but LFAM 5 and GSS 30 IRS-3 are confirmed as radio sources. 

5.2. X-ray Fluxes and Spectral Pcirameters 

The absorbed fiuxcs in column 8 of Table 1 were computed from the event list of each 
source, which contains the time-tagged energy of each detected photon. An auxiliary response 
file (ARF) was generated for each source, which provides the effective area corrected for 
quantum efficiency (cm^ counts photon"^) versus energy (keV) at the source's position on 
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the detector. The absorbed X-ray flux at Earth is, 

1 " E- 

= 7 E ^ - /bkg^, (1) 

where is the energy of each photon, ai{E) is the corresponding effective area, t — 
96439.96 s is the hvetime of the observation, A is the area of the source extraction region, and 
/bkg = 3.605 X 10~^^ ergs cm~^ s~^ pixel"^ is the 0.5-7.0 keV background flux as determined 
from a deep, source-subtracted observation from the Chandra calibration database. 

Traditionally, extraction and model fitting of a background-subtracted source spectrum 
are used to determine spectral parameters such as the absorption column density Nh, mean 
plasma energy kT, and unabsorbed X-ray luminosity Lx- Obviously, an accurate determi- 
nation of unabsorbed Lx requires knowledge of Nh, as discussed by Casanova et al. (1995). 
For some sources in p Oph, reliable estimates of the visual extinction are available from 
previous studies, and we are able to constrain the absorption column density using the re- 
lation Nb/Ay = 1.57 X 10^1 cm-2 mag"^ (Vuong et al. 2003) (§2). In the absence of Ay 
information, the X-ray data were used to estimate Nh. 

Although the traditional spectral fitting approach gives satisfactory results for bright 
sources ( > 100 counts), it does not work well for fainter sources because of issues related 
to binning of the spectrum when only a few counts are available. Since ~60% of the X-ray 
detections in p Oph A have <100 counts, we have developed a non-parametric method for 
estimating spectral parameters that operates on unbinned photon event data, as described in 
more detail in Appendix B. This method is based on a large number of simulated isothermal 
ACIS-I spectra spanning a range of (logiVn, logT) values, utilizing the APEC emission 
model (Smith, Brickhouse, Liedahl, & Raymond 2001), the Wisconsin absorption model 
(Morrison & McCammon 1983), and the MARX ray-tracing Chandra simulation software 
^. The method was used to determine the X-ray temperatures and unabsorbed luminosities 
for all Chandra detections given in Table 4, as well as Nh values for those detections where 
a previous determination of Ay was not available. The X-ray temperature and luminosity 
estimates are based on the assumption that the plasma is isothermal. The inferred values 
thus provide only a rough temperature estimate that can be used to compare a large number 
of sources. In reality, the plasma is likely to be more complex, time variable and involve 
multiple temperature components. 



'5http://space.mit.edu/CXC/MARX/ 
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5.3. X-ray Variability and Light Curves 

To quantify X-ray variability, the event list of each source was used to calculate the KS 
statistic, given by 

KS = V^sup|/,(t)-/o(t)|. (2) 

where n is the number of events, /j(t) is the normalized observed cumulative distribution 
and /o(t) is the normalized model cumulative distribution, assuming a constant flux. If z 
is the value of the KS statistic then the null hypothesis probability P(-2), or probability of 
constant flux, is given by a convergent infinite series (Press et al. 1992): 



nz) = 2Y,{-iy-'e-'^''' (3) 
i=i 

Truncating the series at six terms gives a few representative values: P(0.5) = 0.96, 
P(l.O) = 0.27, P(1.5) = 0.02, P(1.7) = 0.006, and P (2.0) = 6.7 x lO'l 

Sources with z — KS > 1 generally show some light curve variability and sources with 
KS > 2 typically reveal obvious signs of flaring. The source event flies were used to generate 
light curves for all 87 X-ray sources. Figures 2-8 show light curves for seven sources showing 
interesting variability. The time bin sizes in each light curve are varied so as to include 
approximately the same number of counts per bin. 

The top panel in each light curve shows fx- which as already noted is the the back- 
ground subtracted flux at Earth (Sec. 5.2). In the simplest approximation, logT and log A^h 
determine the shape of the X-ray spectrum and the flux conversion factor, the ratio of 
unabsorbcd-to-absorbcd flux. At each time bin, logT and logA^n and their associated ran- 
dom errors are estimated by the non-parametric method described in Appendix B. The two 
lower panels in Figs. 2-8 show logT and unabsorbed luminosity Lx assuming a distance 
d — 165 pc. 

Fig. 2 reveals that the optically bright WTTS DoAr 21 was in the decay phase of a 
large flare throughout the observation, and the spectral changes that occurred during this 
decay are discussed below (Sec. 7.2.3). As apparent in Fig. 2, the decay was punctuated 
by frequent hard-energy spikes, which can be interpreted as reheating events. The peak in 
mean energy often coincides with or immediately precedes the luminosity peak. The time- 
resolved spectral flts suggest that the hard-energy spikes are not consistent with sudden 
increases in column density. As a result, it may necessary to explicitly consider reheating 
when modeling flare decay with loop models. The relatively short rise and decay times 
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suggest moderate-sized magnetic loops. 

5.4. Extragalactic X-ray Sources 

The KS statistic, count rate, and (E) can be used to isolate possible extragalactic X-ray 
sources (Daniel, Linsky & Gagne 2002), as demonstrated in Figure 9a-c. Of the 27 Chandra 
sources with no IR or radio counterparts, 26 show (E) > 3 keV, count rates < 0.001 cts s~^, 
and KS< 1. XSPEC simulations of faint, highly absorbed (25 < Ay < 50) power- law 
(2 < a < 4) sources sources give rise to a similar distribution of (E) and count rate. 
Furthermore, extragalactic X-ray sources are expected to have low KS statistic because they 
generally do not show large-amplitude variations on short time scales (Zamorani et al. 1984; 
Cihegi & Maccacaro 1997). 

Using equations (2) and (3) of Cowie et al. (2002) to estimate the expected number of ex- 
tragalactic sources in a 96 ksec exposure of p Oph, a limiting flux /x = 9 x 10~^^ ergs cm^^ s~^ 
yields approximately 16 extragalactic sources in the 16^9 x 16'9 AClS-1 FOV. We note, 
though, that the predicted extragalactic source count depends on the assumed extragalactic 
source spectrum and the spatially varying absorbing cloud column, both of which are diffi- 
cult to estimate. Nonetheless, it seems hkely that the majority of the 26 Chandra sources 
with no infrared or radio counterparts are extragalactic. However, a few of the 26 uniden- 
tified sources may be new cloud members. In particular, the unidentified Chandra source 
J162627. 4-242418 is a candidate for cloud membership since it is variable in X-rays (KS = 
2.77) on a timescale of less than one day. 

5.5. X-ray vs. Infrared Properties 

Figure 10 shows a, J — K versus K color-magnitude diagram for those IR sources in 
the ACIS-I FOV with measured JHK magnitudes, including a few fainter sources which 
have only lower magnitude limits at J. The faintest K-band source detected by Chandra is 
J162651. 9-243039 (= GY 144) at K ^ 13.46, although its J color is not well-determined. 

In Figure 11 we show a color-color diagram for those IR sources having measured J, H, 
and K magnitudes. As Fig. 11 shows, most of the Chandra detections are confined to a strip 
oriented roughly in the direction of the Ay = 10 reddening vector, consistent with normal 
reddening for cluster members. The colors of the CTTS GY 81 and the WTTSs SKS 1-11 
and GY 84 are consistent with cluster membership but they were not detected with Chandra. 
We note that the source GY 101 has unusual colors for its assigned class III type, but its 



J-H value is only a lower limit. We used J > 17.0 (Table 4) to compute its J-H in Fig. 11, 
but 2MASS gives a more stringent limit J > 18.67. 



5.5.1. Unclassified Infrared Sources 

Chandra detected X-ray emission from more than two dozen IR sources whose SED 
classifications are not yet known (Table 4). Are these unclassified IR sources cloud members, 
such as T Tauri stars? We argue below that most of these X-ray /IR sources are not X-ray 
active AGNs, and conclude that many are hkely to be cloud members. 

To substantiate the above conclusion, we consider the ratio of X-ray to K-band fiuxes. 
To estimate the unabsorbed IR fluxes, we calculated the interstellar extinction for the J, H, 
and K bands using 

2.4A^-^Vm, (4) 



E(J-ir) 

as derived by (Draine 1989). In conjunction with the fact that. 



E{J-K) 



1.1 X 10^^cm"^mag-\ (5) 



(Draine 1989), we derived Ax as a function of known Ay for each source. We found that 
for a typical source {kT 2.0 keV and A^h = 3.27 x 10^^ cm~^), the K band extinction 
(/k//ko — 0.15) is similar to the X-ray extinction {fx/ fxo — 0.25). Therefore, we used 
/k to analyze the relationship between IR and X-ray emission. Figure 12 illustrates that 
log(/x//K) is slightly lower for unclassified X-ray/IR sources (median log/x//K = —3.0) 
than YSOs with measured SEDs (median log/x/Zx = —2.6). This suggests that most of the 
unclassified X-ray/IR sources arc probably not X-ray active AGN because AGN typically 
show higher X-ray-to-optical flux ratios than normal stars (Maccacaro et al. 1988; Stocke et 
al. 1991). 



6. Radio Source Properties 

6.1. Radio Variability 

Each of the VLA detections in Table 3 was observed in two or more ~15 minute scans 
separated in time by anywhere from ~0.5-3 hours. To search for short-term (~hours) vari- 
ability, we compared peak fluxes in the individual scans for the brightest radio detections 
(S/N>10). 
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Definite variability was detected in LFAM 2 (J162622.4-242252) and Oph SI (J162634.2- 

242328). The peak 6-cm flux of LFAM 2 increased by a factor of ~5 during the time interval 
0629 - 0931 UT, indicative of a radio flare. This radio source is associated with the class III 
IR source GSS30 IRS-2 (Barsony et al. 1997). No large amplitude flares were detected from 
Oph SI, but its peak fluxes showed significant scan-to-scan variations of up to ~30% {30a) 
over time intervals of ~3 hours. For example, a 6 cm scan centered on Oph SI in the time 
interval 0629 - 0644 UT gave a peak flux S^^^ = 7.94 ± 0.06 mJy (ilcj) while a second 
scan at the same pointing position from 0921 - 0931 UT gave = 5.86 ± 0.08 mJy. Four 
additional scans with Oph SI positioned 6'A off-axis also confirm a trend of decreasing radio 
flux during the observation. Thus, we find for the first time that radio variability is present 
in this unusual close binary system. We note little or no direct correlation between the VLA 
and Chandra variability of Oph SI. For example the X-ray flux peak from 0820 - 1005 UT 
does not coincide with the peak 6-cm emission interval from 0629 - 0644 UT. The apparent 
decoupling between the X-ray and radio emission is not unique to Oph SI, and has been 
observed in other PMS objects such as the multiple sytem V773 Tau (Sec. 7.3). If the radio 
emission is produced in an extended magnetosphere around the B4 star (Andre et al. 1991) 
and the X-ray flares are produced in a corona around the K-type secondary, then one would 
not necessarily expect to see a correlation between the X-ray and radio emission. 

Given that the WTTS DoAr 21 was in the decay phase of a large X-ray flare, it is 
noteworthy that its 6-cm radio flux appeared rather stable. A 15-minute scan centered at 
0619 UT with DoAr 21 located 2(7 from phase center gave a peak flux of S^^^ = 12.04 ± 
0.05 mJy, and a second 10-minute scan at the same pointing position centered at 0939 UT 
gave ^fP'f = 11.82 ± 0.04 mJy Two additional scans centered at 0542 UT and 1007 UT 
with DoAr 21 located 7'7 from phase center also gave peak fluxes differing by less than 3o". 
Thus, we find no compelling evidence for radio variability in DoAr 21 up to the end of the 
last VLA 6 cm scan at 1013 UT, but no simultaneous radio coverage is available during the 
last ~16 hours of the Chandra observation when several X-ray temperature spikes occurred 
(Fig. 2). Even though no 6 cm variability was detected, DoAr 21 is clearly variable on longer 
timescales. The 6 cm fluxes quoted above are between the peak flare value of 46.1 mJy on 
1983 February 18 and the much lower level of 3.7 mJy measured on 1983 June 18 (Stine et 
al. 1988). 

6.2. Circulcir Polctrization 

Circularly polarized radio emission indicative of nonthermal processes was clearly de- 
tected at 6 cm in three sources. These are the optically visible G-K type star DoAr 21 
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(J162603.0-242336), Oph SI (J162634.2-242328), and the extragalactic source LFAM 21=BZ6 
(J162700. 0-242640). The fractional circular polarization (tTc =V/I) and detection significance 
(S/N ratio) for these sources are given in the notes to Table 3. 

Previous detections of circular polarization in Oph SI at both 6 cm and 2 cm (with an 
upper limit at 20 cm) were reported by Andre et al. (1988). We confirm the presence of 
circular polarization at 6 cm and also provide a new detection at 3.5 cm, but likewise fail to 
detect circular polarization at 20 cm. 

7. Discussion: The Young Stellar Population in p Oph A 

7.1. CIeiss O/I Protostars 

The core-A region contains three well-known protostars, namely the class source LFAM 
5 and the class I sources GSS30-1RS1 and GSS30-1RS3. In many respects, LFAM 5 is the 
most unusual of the three. It is a centimeter radio source that is thought to be driving a 
molecular outfiow (Andre et al. 1990), but has so far not been detected at miUimeter, IR, or 
X-ray wavelengths. Our VLA observation confirms radio emission from LFAM 5 (Table 3) 
at a 6 cm fiux density similar to that reported by LFAM. 

GSS30 is an asymmetric infrared bipolar reflection nebula with at least 3 IR sources 
seen in projection (Grasdalen et al. 1973; Castelaz et al. 1985; Weintraub et al. 1993). The 
class I source GSS30-IRS1 is thought to be the illuminating star (Castelaz et al. 1985). It has 
not been previously detected at centimeter radio wavelengths and we likewise obtain only an 
upper limit of Seem < 0.12 mJy (3o") and an identical Tipper limit at 3.5 cm. Observations at 
2.7 mm have also failed to detect GSS30-IRS1 (Zhang, Wootten, & Ho 1997). There is no 
confirmed detectable outflow in this source. However, line emission in recent VLT IR spectra 
was interpreted as possible evidence for dense gas in an accretion shock (Pontoppidan et al. 
2002). 

The class I source GSS30-IRS3 has K = 12.85 but no emission peak is seen at J or H 
(WcintraTib ct al. 1993), and is thus much fainter than IRSl. A 2.7 mm continuum source 
has been detected toward GSS30-IRS3 at 82.7™,™ = 29 mJy (Zhang, Wootten, & Ho 1997). 
Gentimeter radio emission is also detected from the radio source LFAM 1, which is identifled 
with GSS30-IRS3. We confirm radio emission from LFAM 1 (Table 3) at a 6 cm fiux near 
the value reported by LFAM. 

None of the three class /I sources discussed above was detected in our Chandra observa- 
tion. Upper limits (Table 5) suggest that their luminosities are below logLx ~ 28.5 ergs s^^ 
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However, these upper limits are based on uncertain Ay and kT. In the absence of specific 
information, we have assumed kT ?s 2.5 keV and logA^H ~ 22.5 cm~^, typical of X-ray de- 
tected TTS in core A. An extinction estimate Ay ~ 38 is available for GSS30-IRS1 (Castelaz 
ct al. 1985), and using this value along with kT ^ 2.5 keV and a 7 count Chandra detection 
threshold wc obtain logLx < 28.4 ergs s~^ (0.5 - 7 keV). Assuming a higher temperature 
kT Q keV, as measured for two class I sources in the cluster IC 348 (Preibisch & Zinnecker 
2002), then the upper hmit for GSS30-IRS1 becomes logLx < 28.2 ergs s^^ 

Wc thus conclude that the X-ray luminosities of GSS30-IRS1 (and also very likely 
GSS30-1RS3) are well below the value of Lx ~ 10'^° ergs s~^, typical of class 1 sources 
detected by Chandra such as those in IC 348 (Preibisch & Zinnecker 2002) and p Oph core- 
F (Imanishi, Koyama, & Tsuboi 2001). However, the rac?io luminosities of both LFAM 5 and 
GSS30-IRS3 are comparable to those of CTTS detected in core A (Sec. 7.2.2). Several fac- 
tors may infiuence the X-ray detectabihty of class O/I sources. Most importantly, variability 
seems to be a factor since several of the class I objects detected by Imanishi, Koyama, & 
Tsuboi (2001) in core F were caught during X-ray fiares. Extinction also plays a role, and 
class objects in particular are believed to have very high absorption that efficiently masks 
any soft X-ray emission that is present. Extinction will also make it more difficult to arrive 
at a complete picture of X-ray production mechanisms in protostars. GSS30 IRSl is a case 
in point. The suggestion that an accretion shock is present (Pontoppidan et al. 2002) raises 
the possibility that soft X-rays may be produced in the shock region, similar to what may 
be occurring in the CTTS TW Hya (Kastner et al. 2002). Since shock- induced X-rays are 
expected to have characteristic temperatures below ~1 keV, such emission (if present) would 
be heavily absorbed and difficult to detect in protostars with typical extinctions Ay ~ 40. 



7.2. T Tauri Stars 

7.2.1. X-ray Luminosity Functions 

Using the unabsorbed Lx values of X-ray detections (Table 4), the median X-ray lumi- 
nosity of the 12 detected CTTS is logLx = 29.92 ergs s~^ which is similar to the median 
for 15 detected WTTS of logLx = 29.79 ergs s~^. To make further comparisons between 
the X-ray luminosities of CTTS and WTTS, we take X-ray upper limits of the undetected 
TTS into account (Table 5). Including both detections and non-detections, all four versions 
of the generahzed Wilcoxon test in the ASURV statistical analysis package (Lavalley et al. 
1992) give a probabihty p — 0.53 that the CTTS and WTTS X-ray luminosities are drawn 
from the same distribution. Thus, the Wilcoxon test is inconclusive for this sample. Even 
so, the the Kaplan-Meier estimators for the CTTS and WTTS samples shown in Figure 13 
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appear very similar, agreeing to within the error bars. 

We also find that the spectral hardnesses in CTTS and WTTS are nearly identical. 
Using the mean photon energies (E) in Table 1, the CTTS detections have median (E) = 
2.57 and WTTS detections have median (E) = 2.53. 

On the basis of the shapes of their IR spectral energy distributions, CTTS (class II 
sources) are thought to be surrounded by circumstellar disks but WTTS (class III sources ) 
are not. Given that the X-ray luminosities and spectral hardnesses of the CTTS and WTTS 
in our sample are quite similar, there is no compelling evidence that the presence or absence 
of infrared disks has any significant effect on their X-ray emission. However, it should be 
emphasized that accretion diagnostics (like Ho; or Ca II emission) arc not generally available 
for sources in p Oph cloud core A. Thus, our results cannot be used to determine whether 
accretion itself affects the X-ray properties of this small sample. 

7.2.2. Radio vs. X-ray Luminosities 

Figure 14 shows the 5 GHz radio luminosities (I/6cm) and X-ray luminosities (I/x) for 
those 6 TTS (class II and III) that were simultaneously detected with the VLA and Chandra, 
along with Chandra upper limits for class O/I sources. We have also included Oph SI, whose 
IR classification is uncertain because of the close binary separation. Even though the sample 
is small, we have checked for a correlation between Lgcm and Lx amongst TTS. Excluding 
Oph SI, the generahzed Kendall's tau test in the ASURV software package gives a correlation 
probability p = 0.15 and the Cox proportional hazard model gives p — 0.64. If Oph SI is 
included (on the suspicion that its faint companion Oph SIB is a TTS), then these values 
increase shghtly to p = 0.35 and p — 0.81, respectively. Thus, we find that radio and X-ray 
luminosity are not strongly correlated in this small sample of detected TTS. 

We do note however that in our small sample there is a tendency for WTTS to be more 
luminous radio sources than CTTS, even when comparing objects with similar Lx. In Table 
3, all four WTTS (class HI) have larger Lgcm values than the two CTTS (class II). This may 
refiect different radio emission processes in WTTS and CTTS, with nonthermal emission 
being more common in WTTS. It can also be seen from Figure 14 that Oph SI is clearly 
situated in the region of higher Lgcm occupied by WTTS. 

Two objects of some interest are the class II sources Elias 24 (J162624. 0-241613) and 
GY51 (J162630. 5-242256), both of which were detected simultaneously with Chandra and the 
VLA. As Figure 14 shows, they have similar Lx, but GY 51 is nearly an order of magnitude 
more luminous in the radio. GY 51 shows a weak X-ray flare near 1115 UT on May 16, while 
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Elias 24 experienced a moderately strong impulsive X-ray flare near 2400 UT. The 3.6 cm 
and 6 cm radio fluxes of GY 51 are nearly identical, suggesting a flat (or non-rising) spectral 
energy distribution and possible non-thermal radio emission. 

The faint M2 star Ehas 24 was only detected at 6 cm, so no deflnitive estimate of its 
radio spectral index is available. However, there is no signiflcant 3.6 cm emission in a scan 
with Elias 24 positioned 6.6' from phase center, and we obtain an upper limit Ss.gcm < 0.2 
mJy (5o"). This upper limit is approximately equal to the measured flux at 6 cm (Table 
3), suggesting a non-positive spectral index a < (Sj, oc z/"). This may be an indication of 
nonthermal emission from magnetically-trapped particles. Such emission would be consistent 
with the impulsive X-ray variability, which is an indicator of magnetic activity. However, 
the radio interpretion is not yet clear since Gomez et al. (2003) have identified Elias 24 as 
the possible driving source of a near-IR knot located 144" away. If this identification is 
correct, then the 6 cm radio emission could be thermal emission associated with an outfiow, 
as is already known to occur in other sources associated with Herbig-Haro objects (Bieging, 
Cohen, & Schwartz 1984). Higher sensitivity multi- wavelength observations will be needed 
to determine if the newly-detected radio emission of Elias 24 is thermal or nonthermal. 



7.2.3. Time-Resolved X-ray Spectroscopy of a Large Flare on DoAr 21 

DoAr 21=GSS 23=V2246 Oph (Kl, WTTS) is the brightest X-ray source in the ACIS-I 
FOV. As is evident in Fig. 2, DoAr 21 was in the decay phase of a large flare throughout 
the Chandra observation. DoAr 21 was placed far off-axis to migitate pile-up in the ACIS 13 
CCD. Although the total count rate of DoAr 21 exceeded 1 counts s~^ at the beginning 
of the observation, the source did not pile up signiflcantly, contrary to a previous claim 
by Imanishi et al. (2002). We flnd that the pileup fraction never exceeded a few percent 
because the maximum count rate in the central pixel of the PSF was 0.015 counts s~^ pixel"^. 
Although the count rate decay in the top panel appears to be quite smooth up to a small flare 
around day 0.8, the logT light curve in the second panel is variable throughout. We note 
that the temperature spikes appear to signal reheating events. In many cases, a temperature 
spike immediately precedes a /x maximum. Based on its V-band absorption (Wilking et al. 
2001), we expect DoAr 21 to have a column density A^h ~ 1-0 x 10^^ cm^^. logT and Lx in 
each time bin were estimated using the non-parametric method in Appendix B. 

In order to conflrm the validity of the light curve analysis, the ~81000 counts associated 
with DoAr 21 were divided into 41 time segments each containing approximately 2000 counts. 
For each of the 41 ACIS-1 spectra, background spectra and response flies were generated 
using the CIAO tool PSEXTRACT. Each of these were modeled in XSPEC VI 1 using a 
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single-temperature VAPEC model. We generally found good fits (x^ ^ 1) using an Fe 
abundance of 0.5 solar and a column density iVn ~ 1.0 x lO^^cm"^. Because column-density 
variations have been reported during X-ray flares on YSOs (Tsuboi ct al. 1998), the column 
density, temperature, and normalization were free parameters. After each fit, the emission 
measure and absorption-corrected X-ray luminosity were calculated assuming a distance of 
165 pc (Chini 1981). Indeed, apparent variations of 10-20% in the best-fit column density 
values are seen throughout the observation. The increases in A^h do not generally show 
corresponding increases in hardness ratio or count rate. The increases in best-fit A^h are, 
however, associated with decreases in the best-fit temperature and increases in Lx as a result 
of the increased correction due to absorption. We believe that this result is an artifact of the 
fitting procedure: moderate-resolution, moderate signal-to-noise spectra cannot distinguish 
small increases in column density from small increases in temperature. Since we have no 
reason to suspect large column-density changes on a WTTS like DoAr 21, the column-density 
was fixed at the median best-fit value: A^h = 1-08 x lO^^cm"^ and the spectra were then 
refit. Note that this value of Ah determined from spectral fits differs by only 6% from the 
value determined by our non-parametric analysis (Table 4). 

Assuming V = 13.62, BC = -0.36, d = 165 pc, A^h = 1-08 x lO^^cm'^, and A^hMv = 
1.57 X 10^^ cm~^ mag~^, the bolometric luminosity of DoAr 21 is Lboi ~ 2.3 x 10^^ ergs s~^. 
During the Chandra observation, Lx/L^oi decreases from 4.0 x 10~^ to 2.2 x 10~^. 

The results of the time-resolved spectral analysis arc shown in Fig. 2. There is a steady 
overall decline in emission measure and Lx over the course of the observation and the plasma 
undergoes no fewer than five temperature spikes. If the temperature spikes are ignored, then 
the plasma would appear to be cooling slowly from a giant flare which occurred before the 
start of the Chandra observation. The long decay time and the slow but steady decrease in 
temperature would indicate a large, low-density flaring loop. For example, the two-ribbon 
flare model of Kopp & Poletto (1993) implies a loop radius in the range 1 — 3Rq. 

Previous observations of DoAr 21 with ROSAT and ASCA suggest approximately one 
large flare per day (Imanishi et al. 2002). If the flare decay in Fig. 2 is typical of flares on 
DoAr 21, then the corona is in a state of nearly continous flaring with low- and moderate- 
energy flares outnumbering large flares. 

The frequent increases in kT seen in Fig. 2 suggest that plasma is being intermittently 
reheated by moderate-energy flares. If the reheating events occurred in the same loops as 
the primary decaying flare then a simple cooling-loop model will signiflcantly overestimate 
the cooling time, thereby underestimating the density and overestimating the loop radius. 
Or, the smaller flares may signal reconnection events in separate magnetic loops. 
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Long, uninterrupted monitoring of DoAr 21 or simliar stars may provide a clue to the 
location of these smaller, shorter-duration flares. If short-duration flares during the decay 
phase of long-duration flares are more frequent than short- duration flares during "quiescent" 
periods, then the two types of flares probably occur in the same or adjoining arcades of 
magnetic loops. 

A final question concerns the geometrical relationship between the X-ray and radio 
emitting regions in DoAr 21. As already noted (Sec. 6.1), two separate VLA scans of 
DoAr 21 at 0619 UT and 0939 UT (16 May) gave nearly identical 6 cm fluxes differing by 
<2%, while at the same time its X-ray flux declined by about 15%. This indicates that the 
circularly polarized nonthermal radio emission originating in large magnetospheric structures 
is effectively decoupled from the thermal X-ray emission, which originates much closer to 
the star - presumably in hot coronal plasma. Similar behavior was noted by Feigelson 
et al. (1994) in simultaneous RO SAT and VL^/ observations of the unusual WTTS 
V773 Tau. However, the radio flux of V773 Tau declined while its X-ray emission remained 
steady, whereas the opposite behavior is observed in DoAr 21. The interpretation of the 
behavior of V773 Tau is now complicated by the discovery that it is a quadruple system 
consisting of a K2V -|- K5V spectroscopic binary with a nearby M-dwarf at 0.12" (Ghez et 
al. 1995), and an IR companion 0.21" away (Duchene, Ghez, & McCabe 2001). ROSAT 
could not spatially resolve these four components, raising the obvious question of whether 
the (variable) radio emission and (steady) X-ray emission originated in the same star. Even 
though DoAr 21 is not yet known to be such a closely-spaced multiple system, this possibihty 
should not yet be dismissed given its similarities to V773 Tau. 



7.3. The Unusual Close Binary System Oph SI 

Figure 15 shows the ACIS-I spectrum of the unusual magnetically-active binary system 
Oph SI. As we have noted, this star consists of a B-type primary and a fainter secondary 
at 20 mas separation (Simon et al. 1995). At this close separation, the two components 
cannot be spatially resolved by Chandra or the VLA. The spectrum shows strong absorption 
with little or no detectable emission below ~1 keV. There are no clearly discernible emission 
line features in the spectrum. The light curve shown in Fig. 4 shows Lx variations from 
2 - 4 X 10^0 ergs s"^ with two small flare peaks at 0815 UT and 2453 UT. 

We have analyzed the spectrum using a variety of emission models in XSPEC includ- 
ing single-temperature (IT) and two-temperature (2T) optically thin plasmas (VAPEC), 
bremsstrahlung, power-law (PL), and the Chebyshev polynomial differential emission mea- 
sure (DEM) model C6PVMKL. All models included an absorption component and the spec- 
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trum was rebinned to a minimum of 15 counts per bin prior to fitting in XSPEC. We 

compared spectral fits using spectra extracted with CIAO v. 2.2.1 and with the more re- 
cent CIAO V. 2.3/CALDB 2.0. and CIAO v 3.0.2/CALDB 2.26 releases. Corrections for 
charge transfer inefficiency (CTI) were applied to the spectra extracted with CIAO v. 2.3 
and CIAO v 3.0.2. The best-fit values from spectral fits in Table 6 are from CIAO v. 2.3 
with CTI corrections applied. Generally, the three different versions of CIAO gave similar 
best-fit values but the fit residuals as measured by reduced were 15% - 20% smaller with 
CTI corrections applied. 

Because of the lack of strong emission features, several different models can acceptably 
fit the Oph SI spectrum. Table 6 gives the derived spectral parameters for a power- law 
model as well as IT and 2T VAPEC models. As can be seen, the fit statistics are very 
similar for these three models. Thus, we cannot distinguish between a power-law spectrum 
and optically thin plasma emission at the signal-to-noise ratio of the Oph SI X-ray data. 

The IT VAPEC model imphes a mean plasma temperature kT 2.4 keV. A bremsstrahlung 
model gives a nearly identical temperature kT 2.1 keV and the C6PVMKL model also 
shows a strong emission measure peak at kT ^ 2.0 - 2.4 keV. Considered collectively, these 
models indicate that most of the emission arises from plasma at or slightly above kT 2 
keV (T 2.3 MK), assuming a thermal origin. 

A IT VAPEC model using solar abundances gives a relatively poor fit (x^/dof = 
165.4/135 = 1.23). The fit can be substantially improved by allowing the Fe abundance 
to vary. As Table 6 shows, the inferred Fe abundance is quite low, amounting to only a few 
tenths of the solar value. Little or no further improvement was obtained by allowing the 
abundances of other elements (Ne, Mg, Si, S) to vary. We thus conclude that if the emission 
is thermal then the iron abundance is well below solar. 

All models give an absorbed fiux log/x = —12.52 ergs cm~^ s~^ (0.5 - 7 keV) and an 
absorption column density near log Nh = 22.27 cm~^. The IT and 2T VAPEC models 
(Table 6) give unabsorbed luminosities in the range log Lx = 30.37 - 30.59 ergs s~^ (0.5 - 7 
keV) at d = 165 pc. The above values are in very good agreement with those derived using 
the independent non-parametric approach (Tables 1 and 4; Appendix B). 

The temperature range deduced from VAPEC models of Oph SI is similar to that seen 
in T Tauri stars. The X-ray luminosity of Oph SI is somewhat higher than the median value 
for TTS detected in core A, but still within the observed range for TTS. Thus, most of the 
X-ray emission could be produced by the faint companion star Oph SI B discovered by lunar 
occultation. 
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7.4. Brown Dweirf Candidates 

Chandra detected 4 of 15 known brown dwarf candidates (BDCs) in core A, namely GY 
5 (spectral type M7), GY 31 (M5.5), GY 37 (M6) and GY 59 (M6) (Wilking et al. 1999). As 
a class, these objects are very weak X-ray emitters. The median luminosity of the four BDC 
detections is logLx = 28.3 ergs s~^, which is a factor of ~50 below the median Lx for TTS 
detected in core A. Assuming a detection threshold of 7 counts, the 11 undetected BDCs 
have upper limits logLx < 27.80 ergs s~^. Using the bolometric luminosities of Wilking et 
al. (1999), the four detected BDCs have logLx/I/boi in the range —4.71 to —3.46, where the 
latter value is for the flare source GY 31. Brown dwarfs in other star-forming regions exhibit 
similar X-ray activity levels (Imanishi et al. 2002; Preibisch & Zinnecker 2002; Feigelson et 
al. 2002; Mokler & Stelzer 2002). 

In addition to their faint X-ray emission, the BDCs in core A are also undetected in 
our VLA observations. Typical radio upper limits for X-ray detected BDCs are S'ecm < 0.16 
mJy or Lgcm < 15.7 ergs s~^ Hz~^. 

Apart from their faint Lx, the X-ray properties of detected BDCs appear to be similar 
to TTS. The median value of mean photon energy for the detected BDCs is (E) = 2.3 keV, 
only slightly lower than the median of 2.5 keV for TTS (§7.2.1). Strong variability in the 
form of a flare was detected in one BDC (GY 31), as discussed below. Thus, GY 31 has a 
larger X-ray to K-band flux ratio fx/ fx than most TTS, but the other 3 BDCs which did 
not flare have fx/fK similar to TTS (Table 4). 

Figure 16 shows the flaring light curve of GY 31. The flux declined by a factor of ~ 5 
during the flrst 10 ksec of the observation. The hard-band light curve (not shown) closely 
tracks the broad-band behavior, but there is very little emission in the soft-band below 2 
keV due to strong absorption. 

We attempted to flt the time-averaged spectrum of GY 31 (Figure 17) using optically 
thin plasma models as well as bremsstrahlung and power-law models. These different models 
give nearly identical reduced values and we conclude that there are insufficient counts 
(Ri380 counts) to rehably distinguish between optically thin plasma emission and other al- 
ternatives. All models give a high absorption column density log Nj^ (cm~^) = 22.76 ± 0.08. 
Equating this to a visual extinction (Vuong et al. 2003) yields Ay = 37 mag. Thus, this ob- 
ject is viewed through very heavy absorption, as also determined from infrared observations 
(Wilking et al. 1999). 

If the emission is assumed to originate in an optically thin plasma as is usually the case 
for TTS, then IT VAPEC fits using a subsolar Fe abundance in the range Fe = 0.3 - 0.5 solar 
give kT — 2.3 [1.6 - 2.9] keV. Due to the high absorption, 2T models with a cooler component 
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give no significant improvement over IT models. This time- averaged temperature structure 
is similar to that found in TTS in our sample. 

8. Summary 

We report new results based on simultaneous X-ray and radio continuum observations 
of p Oph cloud core A using Chandra and the VLA. The most important findings of this 
study are the following: 

1. Chandra detected 87 X-ray sources, of which 60 have known infrared or radio counter- 
parts. At least 16 of the 26 unidentified Chandra sources are suspected to be extra- 
galactic. More than one-half of the X-ray detections showed variability, and the X-ray 
sources are typically hard and heavily absorbed ((E) ^ 3 keV, Ay ~ 2 - 56 mag). 

2. Chandra detected 12 of 14 known CTTSs and 15 of 17 WTTSs in p Oph core A. 
Their X-ray luminosities and spectral characteristics are similar, and we thus find no 

compelling evidence that the presence or absence of circumstellar disks (as discerned 
from IR excesses) has any significant effect on their X-ray emission. The effects (if 
any) of accretion on the X-ray emission in this small TTS sample are not yet known 
due to the lack of suitable (e.g. optical) accretion diagnostics. 

3. None of the three class O/I protostellar sources in core A was detected by Chandra 
with upper limits at logLx ^ 28 ergs s~^. A comparison with X-ray detections of class 
I sources in other star-forming regions suggests that such protostars are preferentially 
detected during periods of enhanced X-ray emission (ffares). Our data indicate that in 
the absence of flares. X-ray luminosities can be one to two orders of magnitude below 
the typical value logLx ~ 30 ergs s~^ reported for flaring class O/I sources. 

4. Chandra detected 4 of 15 known brown dwarf candidates in core A, including the flaring 
source GY 31. Their X-ray emission is faint with a median logLx = 28.3 ergs s~^ that 
is a factor of ~50 lower than T Tauri stars in core A. The mean X-ray energies and 
variability of BDCs appear to be similar to TTS. None of the X-ray emitting BDCs 
was detected with the VLA. 

5. The WTTS DoAr 21 was in the decay phase of a large X-ray flare during the Chandra 
observation, but no radio variablity was detected. Time-resolved X-ray spectroscopy 
reveals multiple secondary flares during the decay. These secondary flares are associ- 
ated with temperature increases, suggesting that the plasma is being reheated if the 
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secondary events occur in the primary flaring structure. Such reheating could con- 
tribute to the long decay times seen in some large YSO X-ray flares. 

6. The X-ray hght curve of the binary Oph SI shows hot, moderate-amphtude, short- 
duration flares typically seen on TTSs. Much of the X-ray emission may be produced 
around the K-type secondary located 20 mas from the B4 primary. It remains to be 
determined whether the variable radio emission originates in a magnetosphere around 
the B4 star or is instead due to the nearby late-type companion. 

7. Multifrequency VLA observations detected 31 radio sources in p Oph core A, of which 
16 are confirmations of previous radio detections. New radio detections are reported 
for the emission-hne star Elias 24 and the optically invisible IR source WLY 2-11, as 
well as the first detection of circular polarization in the radio-bright WTTS DoAr 21. 
There is no significant correlation between the X-ray and radio luminosities in the 
small sample of T Tauri stars detected simultaneously with Chandra and the VLA. 
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A. Electronic Database of Objects in p Oph Core A 

The p Oph A region has been studied extensively in most regions of the spectrum and 
a large amount of observational data exist. In order to make comparisons between the new 
Chandra and VLA data and existing data, we used IDL to construct a database of the 
cataloged infrared (Barsony et al. 1997; Bontemps et al. 2001; Allen et al. 2002, 2MASS^), 
X-ray (Casanova et al. 1995; Grosso et al. 2000) and radio sources (Stine et al. 1988; Leous et 
al. 1991) in the Chandra FOV. We matched sources which fell within the combined positional 
uncertainty published for each source. We visually checked every match using the ds9 tool 
in CIAO and against the alternative associations determined by the author of each catalog. 
Similarly, the measured photometry of each source was also checked for consistency. We 
omitted two sources from the 2MASS catalog as they had no counterparts and appeared to 
be reflection nebulae in the 2MASS images. 



'^2MASS 2002, 2d Incremental Release, Point Source Catalog. 
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We used the findings of Simon et al. (1995), who performed a search for binaries by lunar 
occultation. Their study included the portion of p Oph discussed in this paper. Three binary 
systems, Oph SI, DoAr 24E, and Elias 30, are in the Chandra FOV. In our database, all 
information previously associated with Oph SI is assigned to Oph SI A. We assign the same 
V-magnitude absorption {Am) and column density (A^h) to both Oph SI A and Oph SIB. We 
associate all Chandra X-ray emission with Oph SI B (see §7.3). Because of the close 20 mas 
separation of Oph SI A and Oph SI B, some caution is needed when attributing emission to 
either of the two components. Any further determinations of the origin of cataloged emission 
for Oph SI are at the discretion of the user of the data base. We also associate the Chandra 
X-ray emission in the region of DoAr 24E with the binary companion, DoAr 24E B. We use 
the IR photometry determined by Allen et al. (2002) for this source. All other cataloged 
information is assigned to DoAr 24E A in this database. The user should be aware that 
the IR photometry measurements for DoAr 24E A are likely combined in the IR catalogs. 
Elias 30 is separated in Barsony et al. (1997), and we use their catalog for the IR photometry. 
We use the K-band magnitudes derived by Simon et al. (1995) for each of these six sources. 

A total of 318 previously cataloged objects exist in the Chandra ACIS-I FOV. In our 
observations we add 27 new objects detected with Chandra^ one of which (J162607. 3-242530) 
was simultaneously detected with the VLA (see §4.2). We matched these 345 sources with 
known brown dwarf (BD) candidates (Neuhauser et al. 1999; Wilking et al. 1999), spectral 
types (Martin et al. 1998; Wilking et al. 1999, 2001), and classification designations based 
on studies of spectral energy distributions (SEDs), using classifications in order of preference 
from Wilking et al. (2001); Luhman & Rieke (1999); Andre & Montmerle (1994); Greene et 
al. (1994); Wilking et al. (1989). 

The Stine et al. (1988); Leous et al. (1991); Casanova et al. (1995); Strom et al. (1995); 
Barsony et al. (1997); Martin et al. (1998); Wilking et al. (1999); Neuhauser et al. (1999); 
Luhman & Rieke (1999); Grosso et al. (2000); Bontemps et al. (2001); Wilking et al. (2001); 
Allen et al. (2002, 2MASS) catalogs were used to construct a database of 345 sources in the 
ChandraFOV. The IR photometry was preferenced using Barsony et al. (1997) measurements 
first, then Allen et al. (2002), then 2MASS, and Bontemps et al. (2001) last. For sources 
with cataloged JHK photometry we used, 

/, = F,,^^^AA X 10-°-^^ (Al) 

to determine the absorbed flux in each band. Prom these values, and the associated absorbed 
X-ray flux fx (Table 1), we also calculated the ratio of fx/ fx for each band. 

The Ay for each source was determined by Wilking et al. (2001), Wilking et al. (1999), 
and Grosso et al. (2000), in order of preference. Ay determinations based only on JHK 
photometry (e.g., Bontemps et al. 2001) were not used except in estimating upper limits 
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(§7.1) because these methods assume a dereddened color to extract the color excess and Ay. 
This can lead to substantial uncertainties in Ay- For example, Bontemps et al. (2001) used 
(J — H)q = 0.85, close to the locus of T Tauri stars of Meyer, Calvet, & Hillenbrand (1997). 
We note though that for WTTS and CTTS, (J - H) is in the range 0.1-1.1. A 0™5 error 
in (J — H)o corresponds to a substantial error of 4™5 in Ay. To determine the A'h to each 
source with a measured Ay we used, 

Nh^Ay- 1.57 X 10^^ cm-^ (A2) 

(Vuong et al. 2003). Brown dwarf candidates were gleaned from Neuhauscr et al. (1999) and 
Wilking et al. (1999). The database of information on the 345 YSOs in the Chandra FOV 
is publicly available^ in a single IDL save file. Other supplemental materials are available at 
this URL. 



B. A Non-pcirametric Method for Estimating Spectral Pcirameters 

In the traditional method for estimating spectral parameters from X-ray CCD data, the 
list of photon energies is binned to create an X-ray spectrum. The spectrum is fit by forward 
folding a model spectrum through an appropriate set of spectral and effective-area response 
files. Model parameters are usually derived by minimizing some statistic like or the Cash C 
statistic. One can compute the statistic over a grid of free parameters to estimate confidence 
limits for each parameter. This method of estimating spectral parameters is computationally 
efficient and often quite robust when each spectral bin contains a signicantly large number 
of counts. However, as the number of counts per bin becomes small, the fit results become 
uncertain because (a) the statistic is not a reliable goodness-of-fit statistic for small N, and 
(b) the binning process is arbitrary: there is no established method for choosing bin centers, 
bin widths and minimum counts per bin. 

Some methods have been devised to address the first problem related to small N statis- 
tics. For example, Gehrels (1986) derived a method for estimating upper and lower limits 
for small N. Also, Nousek & Shue (1989) discussed various minimization techniques using 
and the C statistic. Mighell (1999) suggest a more robust statistic for small N: Xj- However, 
these methods do not address the second problem of how to bin the event data. 

The above problems - particularly the second - can be avoided by operating on unbinned 
photon event lists, rather than binning the data according to photon energy. Using this 



^ftp: / / astro.wcupa.edu/ pub/mgagne/roph 
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strategy, one can construct empirical (or cumulative) distribution functions (EDF) and then 
take advantage of a number of non-parametric goodness-of-fit statistics that use EDFs (e.g. 
Stephens (1974); Babu & Feigelson (1996)). We briefly describe below a procedure that 
we have developed for estimating spectral parameters (A'h, logT) using unbinned photon 
data. These parameters are then used to estimate the unabsorbcd-to-absorbed flux ratio 
from which we derive the unabsorbed X-ray luminosity Lx from fx and distance. We used 
this procedure with unbinned photon event lists for sources in p Oph A to derive the spectral 
parameters in Table 4 and Figs. 2-8, along with their la uncertainties. As noted in Table 4, 
these errors do not include systematic effects such as time variability, calibration uncertainties 
in the response files used in the MARX simulations, and uncertainties associated with the 
choice of spectral models. For this paper, we used the Wisconsin absorption model (Morrison 
& McCammon 1983) and a single-temperature VAPEC emission model (Smith, Brickhouse, 
Liedahl, & Raymond 2001) with solar abundances except for Fe set at 0.5 times solar. 
If, for example, the X-ray emission were time-variable or involved multiple temperature 
components or had different elemental abundances, then logT would lie outside the stated 
range. Distance uncertainties would lead to much larger errors in Lx- 

XSPEC was used to generate photon spectra over a grid of column densities (20.30 < 
logiVn < 23.30 in increments of 0.04) and plasma temperatures (6.30 < logT < 8.35 in 
increments of 0.05). These spectra were then input to the MARX ray-tracing software to 
generate simulated event files for 75 sources distributed over the ACIS-I detector (similar to 
the situation encountered in p Oph). We collected 16,000-40,000 counts for each (logA^n, 
log T) pair, from which we extracted event lists for thousands of simulated sources of varying 
brightness in the range 10-1000 counts. We then compared the EDFs of these low-count 
sources with bright samples containing 8192 counts. These comparisons were made using 
both the Cramer-von Mises (CvM) statistic and the two-sample Kolmogorov-Smirnov (KS) 
statistic. The choice of 8192 counts for the bright samples was somewhat arbitrary, the 
main consideration being that the bright source samples contain enough counts to totally 
dominate the noise in simulated spectra. 

In order to test the robustness of the method, we tested the method's ability to recover 
the input model from a fake event list. For example. Figure 18 shows confidence contours for 
a fake source with input logA^H = 22.18 and logT = 7.40. The cross indicates the location 
of the input values. In one realization (upper panel of Fig. 18), 25 source counts (and the 
appropriate number of background counts) are extracted from a larger list corresponding to 
a bright source and its empirical distribution function (EDF) is constructed. The 25-count 
EDF is then compared to the bright-source EDF for all (logA^H, log^) pairs by computing 
the CvM statistic. The dot indicates the logiVn, logT with the lowest CvM statistic CvMmin 
for this realization. The contours indicate CvM-CvM^in = ACvM = 0.174, 0.284 and 0.347, 
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corresponding respectively to 68%, 90% and 95% confidence contours for the CvM statistic 
(Stephens 1974). The bottom panel of Fig. 18 is a similar plot for a 60-count source. Visual 
inspection of these plots and many like them show that the best-fit parameters (dot) are 
generally close to the input values (cross) and lie within the 68% contours approximately 
68% of the time. 

To quantify the uncertainty in the best-fit parameters using this method, 200 simulations 
were realized for each (logT, logA^n) pair. The standard deviation of the 200 output best-fit 
logT and logA^n values was computed for each input (logT, logA^n)- This was done for 10, 
25, 60, 160, and 400-count sources with background addition. This was also done for the two- 
sample KS statistic. Some results of this analysis for the CvM statistic are shown in Table 7. 
The first three columns show the input model: logiVn, logT and FCF, the unabsorbed-to- 
absorbed flux ratio. The next column is source counts. The next four columns show the 
output results: mean and standard deviation of the best-flt log A^h and log T values. These 
standard deviations are shown in Fig. 18 as dashed boxes around the best-fit value for that 
realization. The last two columns in Table 7 show the mean and 68%-confidence upper and 
lower bounds on the FCF used to calculate Lx. 

Table 7 allows us to make the following observations: for moderate column densities 
(logA^H = 21.50 and 22.18) and high temperatures (logT = 7.4 and 7.8), logiVn, logT and 
FCF can be estimated with some precision with as few as 25 counts. If logiVn is very high or 
logT is low, then the parameters cannot be reliably estimated, even with 400 counts. That 
is, photons from a luminous, cool, highly absorbed source will be difficult to distinguish from 
a less luminous, hotter, less absorbed source. Hence the large uncertainities in FCF for large 
A^H and low T. 

The ability to estimate spectral parameters can be further improved by constraining 
logA^H or logT. A similar analysis in which logA^n is known ±10% is shown in Table 8. It 
is clear that if Ay data can be used to constrain A^h, then logT, FCF and hence Lx can be 
reliably estimated in most cases with as few as 25 counts. Finally, we note that although 
both the CvM and KS statistic worked well, the ACvM confidence intervals of (Stephens 
1974) were more accurate than those for the KS statistic for low-count sources. With more 
than 100 counts, both statistics found very similar errors. Thus, in estimating logT and 
Lx in Table 4 and in Figs. 2-8, we used the CvM statistic and pubhshed Ay data when 
available. 

To estimate spectral parameters from real data, photon-event (evt2.fits) and ancillary- 
response (.arf) files are extracted for the source and the source EDF is computed. Based on 
the exposure time and the size of the source region, the number of 0.5-7.0 keV background 
counts in the source region is estimated. The source EDF is then compared to the full set of 
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simulated event lists. Each (logA^n, logT) pair on our VAPEC grid has a simulated event 
file. Based on the number of simulated events, the number of real source events, and the 
estimated number of background counts, the appropriate number of background events are 
randomly drawn from a large background event list and added to the simulated events. The 
simulated background-added EDF is computed and compared to the real source EDF. This 
way, a CvM (or KS) statistic is computed at every point on the (logiVn, logT) grid. A 
full grid search requires approximately 7 seconds on a dual-processor 2.8-GHz Xeon machine 
running IDL version 6.0 under Linux kernel 2.4. 

When log A^H is estimated from Ay^ then the grid search is restricted to values of log A^h 
within 10% of the estimated value. This reduces the error in logT, FCF, and Lx. The best- 
fit values in Table 4 represent the grid point with the lowest CvM statistic. When generating 
hght curves hke those in Figs. 2-8, the procedure is performed at each time step. 

As a reality check, we compared the fit results and confidence contours with those from 
conventional XSPEC fitting. For sources with N > 100 counts, the results were nearly iden- 
tical, provided our VAPEC grid had sufficient resolution. We found that the confidence 
contours for simulated sources with 25 - 100 counts were somewhat smaller. This demon- 
strates that the spectral information content of the unbinned EDF is comparable to or greater 
than that of binned spectra. 

In summary, the advantage of fitting binned background-subtracted spectra in XSPEC 
or Sherpa is that the minimization process is more efficient than deriving confidence contours 
from XSPEC and MARX simulations. The non-parametric method we describe provides 
spectral information for faint sources that may be more statistically valid because the photon 
energies are not binned. 
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Table 1. Chandra X-ray sources in p Oph Core A. 



\JiLJO 




R.A. 


ueci. 


Counts 




ivo 


{£,) 


log /x 




Designation 


(J2000) 


\rt ZUUU j 








yKfi V J 


(org cm s ^ ) 


Designation 






(2) 










I7\ 
\l) 






J162555. 3-242536 


16 


25 


55.40 


-24 25 36.5 


16.7 


6.8 


0.65 


3.56 


-14.24 




J162603. 0-242336'' 


16 


26 


03.03 


-24 23 36.1 


80776.5 


285.4 


16.81 


2.41 


-10.98 


DoAr 21 


J162605. 7-242755 


16 


26 


05.72 


-24 27 56.0 


19.9 


6.3 


0.50 


3.21 


-14.48 




J162607. 1-242724 


16 


26 


07.13 


-24 27 24.3 


731.6 


28.3 


2.51 


3.12 


-12.87 


BKLT J162607-242725 


J162607.3-242530* 


16 


26 


07.36 


-24 25 30.9 


24.9 


6.8 


0.59 


3.97 


-14.14 




J162607.6-242741 


16 


26 


07.63 


-24 27 41.4 


880.7 


30.9 


3.89 


2.88 


-12.85 


BKLT J162607-242742 


J162610. 3-242054 


16 


26 


10.35 


-24 20 55.0 


184.7 


14.8 


3.66 


3.38 


-13.23 


GSS 26 


J162612. 5-242848 


16 


26 


12.51 


-24 28 48.5 


76.5 


10.3 


0.67 


3.60 


-13.73 




J162613. 3-242823 


16 


26 


13.32 


-24 28 23.3 


15.8 


5.7 


0.55 


3.44 


-14.42 




J162615. 7-242749 


16 


26 


15.75 


-24 27 49.4 


58.1 


9.0 


0.59 


3.70 


-13.88 




J162615.8- 241922 


16 


26 


15.83 


-24 19 22.2 


165.7 


14.0 


1.34 


2.42 


-13.71 


SKS 1-7 


J162616. 3-241841 


16 


26 


16.35 


-24 18 41.9 


12.4 


4.7 


0.84 


4.25 


-14.47 




J162616. 8-242223 


16 


26 


16.87 


-24 22 23.0 


2078.2 


46.7 


7.74 


2.47 


-12.59 


Oph S28 


J162617.0-242021 


16 


26 


17.08 


-24 20 21.5 


3772.8 


62.5 


4.23 


1.61 


-12.58 


DoAr 24 


J162617.1-241240 


16 


26 


17.15 


-24 12 40.4 


126.2 


13.5 


1.46 


2.02 


-13.89 


BKLT J162617-241241 


J162617.4^241809 


16 


26 


17.49 


-24 18 09.0 


39.7 


7.6 


0.93 


3.73 


-14.08 




J162617.9- 241518 


16 


26 


17.97 


-24 15 18.6 


41.1 


8.1 


0.82 


3.21 


-14.12 




J162618.8- 242819 


16 


26 


18.88 


-24 28 19.8 


381.7 


20.7 


1.31 


2.54 


-13.29 


VSSG 1 


J162621.5-242600 


16 


26 


21.54 


-24 26 01.0 


10.8 


4.7 


1.24 


2.15 


-15.09 


GY 5 


J162622. 2-242447 


16 


26 


22.23 


-24 24 48.0 


8.8 


4.3 


0.53 


2.59 


-15.01 


SKS 3-12 


J162622.4-242252'' 


16 


26 


22.41 


-24 22 52.8 


3493.1 


60.1 


2.50 


3.04 


-12.25 


GSS .30#2 


J162623.3- 242059 


16 


26 


23.39 


-24 20 59.6 


992.7 


32.6 


3.04 


2.21 


-12.99 


DoAr 24E A 


J162623.6-242439 


16 


26 


23.61 


-24 24 39.6 


445.6 


22.2 


3.52 


4.44 


-12.90 


GY 21 


J162624.0-242448 


16 


26 


24.07 


-24 24 48.1 


2625.2 


52.3 


24.17 


3.45 


-12.24 


Oph S2 


J162624.0- 241613* 


16 


26 


24.09 


-24 16 13.3 


1923.8 


45.0 


3.04 


2.60 


-12.56 


Elias 24 


J162624.1-241518 




26 


24.17 


-24 15 18.7 


20.4 


7.1 


0.88 


3.21 


-14.25 




J162625. 2-242323 


16 


26 


25.24 


-24 23 24.0 


378.0 


20.5 


4.71 


3.52 


-13.14 


GY 31 


J162625. 2-242444 


16 


26 


25.30 


-24 24 44.7 


21.0 


5.8 


1.03 


2.12 


-14.65 


GY 29 


J162625. 7-241427 


16 


26 


25.72 


-24 14 27.3 


105.7 


12.1 


3.16 


1.82 


-14.10 


BKLT J162625-241430 


J162627.4- 242418 


16 


26 


27.48 


-24 24 18.1 


24.9 


6.2 


2.77 


5.18 


-14.02 




J162627.8-242642 


16 


26 


27.83 


-24 26 42.6 


10.5 


4.4 


0.83 


2.00 


-15.00 


GY 37 


J162627.8-242359^ 


16 


26 


27.88 


-24 23 59.1 


12.6 


4.7 


1.91 


5.24 


-14.29 


LFAM 6<= 


J162628. 5-241540 


16 


26 


28.51 


-24 15 40.4 


17.2 


5.4 


0.53 


2.62 


-14.62 


BKLT ,1162628-241543 


J162629.7-241905 


16 


26 


29.70 


-24 19 05.3 


316.8 


18.9 


1.02 


2.87 


-13.33 


SKS 1-19 


J162630.5-242256* 


16 


26 


30.55 


-24 22 56.8 


449.4 


22.3 


3.72 


3.07 


-13.13 


GY 51 


J162630.6- 242023 


16 


26 


30.62 


-24 20 23.3 


22.3 


5.9 


0.68 


4.30 


-14.23 




J162630.8- 242108 


16 


26 


30.84 


-24 21 08.9 


10.0 


4.6 


0.94 


3.65 


-14.58 




J162630. 9-243106 


16 


26 


30.91 


-24 31 06.6 


15.6 


5.9 


0.82 


2.87 


-14.46 


GY 54 


J162631.3-241833^ 


16 


26 


31.32 


-24 18 33.3 


11.5 


4.8 


0.42 


3.59 


-14.48 


LFAM 10 


J162631.3-242530 


16 


26 


31.37 


-24 25 30.3 


36.0 


7.1 


1.26 


2.45 


-14.40 


GY 59 


J162631.5- 242318 


16 


26 


31.52 


-24 23 18.4 


21.6 


5.8 


0.54 


3.97 


-14.33 




J162633.2- 241951 


16 


26 


33.25 


-24 19 51.3 


13.4 


4.8 


0.77 


3.62 


-14.48 




J162634. 1-242328 


16 


26 


34.19 


-24 23 28.2 


2583.1 


51.9 


1.61 


2.52 


-12.51 


Oph SI B 


J162634.5-242228 


16 


26 


34.51 


-24 22 28.7 


39.4 


7.4 


0.96 


3.90 


-14.08 





Table 1 — Continued 



GDS R.A. Dccl. Counts a KS {E) log/x IR/Radio'' 

Designation (J2000) (J2000) (cts) (cts) (keV) (erg cm-^g-i) Designation 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



J162636.8- 241933 


16 


26 36.81 


-24 


19 


33.1 


14.0 


5.0 


0.62 


3.62 


-14.60 




J162636.8- 241551 


16 


26 36.83 


-24 


15 


51.2 


106.2 


11.9 


0.76 


3.49 


-13.64 


BKLT J162636-241554 


J162636.8-241900 


16 


26 36.86 


-24 


19 


00.1 


8.2 


4.1 


0.42 


3.00 


-14.96 


BKLT J162636-241902 


J162637.1-241601 


16 


26 37.12 


-24 


16 


01.1 


64.0 


9.6 


0.76 


1.80 


-14.26 


BKLT J162637-241602 


J162637.9-241942 


16 


26 37.90 


-24 


19 


42.7 


9.4 


4.3 


0.71 


3.45 


-14.83 




J162640. 4-242714 


16 


26 40.48 


-24 


27 


14.4 


162.4 


13.8 


5.02 


3.69 


-13.48 


GY 91 


J162641. 8-241603 


16 


26 41.88 


-24 


16 


03.3 


18.3 


6.0 


0.94 


3.54 


-14.35 




J162642.3- 242625'" 


16 


26 42.39 


-24 


26 


25.9 


45.4 


7.9 


0.75 


3.18 


-14.13 


GY 101 


J162642.5- 242631 


16 


26 42.58 


-24 


26 


31.7 


117.2 


11.9 


0.75 


3.49 


-13.66 


GY 103 


J162642.8-242029 


16 


26 42.89 


-24 


20 


29.9 


1225.2 


36.1 


2.72 


2.54 


-12.80 


GY 110 


J162642.9-242259 


16 


26 42.92 


-24 


22 


59.0 


11.4 


4.6 


1.28 


2.71 


-14.89 


GY 109 


J162643. 4-242434 


16 


26 43.41 


-24 


24 


34.9 


19.6 


5.6 


1.05 


3.88 


-14.42 




J162643.8-241632'' 


16 


26 43.80 


-24 


16 


32.9 


565.7 


25.0 


1.08 


2.35 


-13.18 


VSSG 11 


J162643.9-242717 


16 


26 43.92 


-24 


27 


17.8 


11.9 


4.7 


0.65 


4.52 


-14.35 




J162645.0- 242307 


16 


26 45.06 


-24 


23 


07.7 


71.4 


9.5 


0.80 


3.42 


-13.90 


GY 116 


J162646.8- 241907 


16 


26 46.81 


-24 


19 


07.9 


66.4 


9.3 


1.58 


2.30 


-14.12 


BKLT J162646-241910 


J162648. 3-242834 


16 


26 48.33 


-24 


28 


34.2 


31.5 


7.1 


1.47 


3.22 


-14.16 


BKLT J 162648-242836 


J162648. 6-242839 


16 


26 48.61 


-24 


28 


39.8 


24.5 


6.5 


0.47 


3.51 


-14.29 


GY 128 


J162649. 1-242143 


16 


26 49.18 


-24 


21 


43.8 


18.9 


5.6 


0.59 


3.71 


-14.33 




J162649.2- 242002 


16 


26 49.27 


-24 


20 


02.7 


1348.3 


37.8 


4.42 


2.53 


-12.78 


GY 135 


J162651.9-243039 


16 


26 51.93 


-24 


30 


39.6 


12.6 


5.7 


0.63 


3.23 


-14.57 


GY 144 


J162652.4^243133 


16 


26 52.44 


-24 


31 


33.0 


30.4 


8.0 


0.63 


4.12 


-13.92 




J162653.4-243236 


16 


26 53.49 


-24 


32 


36.7 


31.9 


7.5 


1.85 


3.10 


-14.19 


GY 146 






9ft f\A f\'i 


—24 


zo 


ZO-O 


1 1 fi 












J162654.0-242840 


16 


26 54.10 


-24 


28 


40.4 


11.0 


5.0 


0.50 


3.93 


-14.58 




J162654.4-241554 


16 


26 54.47 


-24 


15 


54.5 


16.9 


6.1 


0.89 


2.63 


-14.56 




J162654.4- 242620 


16 


26 54.48 


-24 


26 


20.6 


428.2 


21.8 


1.79 


2.64 


-13.25 


GY 153 


J162655.0- 242229 


16 


26 55.00 


-24 


22 


29.6 


285.2 


18.0 


1.69 


2.79 


-13.40 


GY 156 


J162658.3-242130 


16 


26 58.38 


-24 


21 


30.2 


41.7 


7.8 


1.11 


3.12 


-14.11 


GY 171 


J162658.6- 241834 


16 


26 58.65 


-24 


18 


34.5 


18.7 


5.7 


0.90 


2.46 


-14.65 


BKLT J162658-241836 


J162658. 9-242215 


16 


26 58.90 


-24 


22 


15.2 


11.5 


4.8 


0.57 


3.34 


-14.70 




J162700. 0-242640'' 


16 


27 00.07 


-24 


26 


40.4 


175.5 


14.5 


0.94 


4.16 


-13.32 


LFAM 21 


J162700. 5-241624 


16 


27 00.53 


-24 


16 


24.9 


58.7 


9.0 


2.41 


1.94 


-14.30 


BKLT J 162700-241627 


J162701.2- 242913 


16 


27 01.28 


-24 


29 


13.6 


12.9 


7.3 


0.72 


3.51 


-14.48 




J162703.0- 242614 


16 


27 03.00 


-24 


26 


14.4 


44.1 


8.1 


2.36 


3.81 


-13.99 


GY 188 


J162704. 1-242828 


16 


27 04.11 


-24 


28 


28.2 


39.8 


8.1 


1.80 


3.28 


-14.11 


GY 192 


J162704.6-242715 


16 


27 04.62 


-24 


27 


15.4 


4614.4 


69.0 


23.33 


3.37 


-12.04 


GY 195 


J162705. 2-242007 


16 


27 05.28 


-24 


20 


07.4 


201.6 


15.5 


1.81 


2.61 


-13.53 


GY 199 


J162706.0-242618 


16 


27 06.01 


-24 


26 


18.5 


228.6 


16.4 


5.82 


3.17 


-13.37 


GY 203 


J162710.3-241912 


16 


27 10.32 


-24 


19 


12.6 


911.4 


31.5 


1.20 


1.59 


-13.22 


Elias 30 A 


J162710.3- 241918 


16 


27 10.35 


-24 


19 


18.5 


110.7 


11.7 


0.98 


1.67 


-14.11 


Elias 30 B 


J162713.7-241816 


16 


27 13.73 


-24 


18 


16.7 


42.3 


8.4 


1.52 


2.42 


-14.31 


VSSG 24 


J162722.9- 241800 


16 


27 22.93 


-24 


18 


00.1 


486.2 


23.6 


2.08 


2.65 


-13.05 


VSSG 22 
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Note. — Col. (1) J2000 lAU designation, cols. (2) - (3) Chandra X-ray position with astrometric correction applied, col. (4) 
net counts (background subtracted), col. (5) Icr error in net counts (Gehrels 1986), col. (6) Kolmogorov-Smirnov (KS) variability 
statistic (variable sources have KS > 1), col. (7) mean photon energy from event file, col. (8) calculated absorbed flux in 0.5 - 7 
keV band (§ 3.1), col. (9) probable counterpart. 

"Simultaneously detected with the VLA. See Table 3. 

'"DoAr=DolidZe & Arakelian (1959); GSS=Grasdalen et al. (1973); VSSG=Vrba et al. (1975); Elias=Elias (1978); LFAM=Leous 
et al. (1991); GY=Greene & Young (1992); SKS=Strom et al. (1995); BKLT=Barsony et al. (1997). 

'^The identification as LFAM 6 is uncertain due to an off'set of R! 2.4" between the X-ray coordinates and precessed LFAM 

coordinates. 
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Table 2. Summary of VLA Observations of p Oph Core A 



Frequency 


Time Range^' 


Pointings 


On-sourcc Time'' 


Flux cal./</i-cal'^ 


Beam FWHM<J 


(GHz) 


(lAT) 




(min) 


(J2000) 


(") 


1.42 


07:41 - 09:16 


1 


20 


0137+331 / 1626-298 


33.0 X 15.9 


4.86 


04:25 - 11:38 




245 


0137+331 / 1626-298 


9.5 X 5.0 


8.46 


05:13 - 08:58 


4 


52 


0137+331 / 1626-298 


5.7 X 2.6 



''Observed on 16 May 2000. 

''Total on-source time was composed of 2 separate scans at 1.42 GHz, 19 scans at 4.86 GHz, and 5 scans 
at 8.46 GHz. 

■^0137 + 331 = 3C48. At low elevations when 1626 — 298 was not visible, phase calibration was achieved 
using 1554 - 270 (4.86 and 8.46 GHz) and 1700 - 261 (4.86 GHz). 

''Typical synthesized beam in cleaned maps; VLA in C-configuration with 26 operational antennas (no 
subarrays) . 

''One pointing centered near Oph SI A and remaining pointings offset by ~one-half primary beam width 
(4.5') in a grid around Oph SI A. 
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Table 3. VLA 6 cm Radio Sources in p Oph Core A ^ 



GDS 


R.A. 
(J2000) 


Decl. 




^ ( Ttpn 

6 cm 

(/«Jy) 


(total) 
6 cm 

(/^Jy) 


RMS Noise 
(/zJy/beam) 


Ident.^ 


Class^ 


J iDZ00D.l-Z4oU14 


1 P. 


'^R in 


94 

Z-i 




1 4 


9 
Z 


4oD 


040 


00 


WLY 2-11 




J iDZDUo.U-Z4ZOO0 


1 R 


9R D"? flR 
ZD UO. UO 


94 
— Z^ 




"iR 
00 


D 


iiyyo 


1 9/1 0/1 
iZ4U4 


4o 


SFAM 7 


777 
111 




1 fi 


9fi 07 Ai^i 


—94 

— Z^i 


Zu 


■^1 

Ol 


Q 


4iU 


^Pi7 
DO 1 


40 








J 10ZOlU.O-Z4zoOo 


1 fi 


9fi in i^ii 


—94 
— z** 


98 
zo 


oo 


Q 
O 




91 1 
Zll 


Q7 
O ( 








71(^0(^11 n O/IOOflQ 
J IDZDi 1 .U-Z4zyUo 


1 fi 
lu 


9R 1 1 (14 
Zu 11. 


94 
— Z^ 


9Q 

zy 


ns 

UO 


z 


A Pi9 
40Z 


/I QO 
4oU 


o4 








J iDZDiO.O-Z4o4Zo 


1 R 


9R 1 t;'? 
ZO 10. Oo 


94 
— Z^ 




98 
Zo 


D 


4uyu 


7/1 7Pi 
/4 / 


79 
1 Z 








TI^^O^^OI 7 0/lOOPil 
J iOZOZi. /-Z4ZZ0i 




9fi 91 7*? 
ZD Zl . / o 


94 
— Z^ 


99 
zz 


t^l 
01 


4 


zoo 


QpiPi 
tjOO 


4o 


LFAM 1 


T 
1 


J iDZDZZ.4-Z4ZZoZ 


1 R 


9R 99 "^Q 
ZD zz.oy 


94 

— Z^ 


99 
ZZ 


00 


A 

4 


9AAI^/'ir^ 

ZDDO^^V J 


9A9n 

ZDZU 


/ 


LFAM 2 


777 
111 


J iDZDZ4.U-Z4iDio 


1 R 


9R 9zL ni 

ZD Z^.Dl 


94 

Z4: 


1 R 

ID 


1 9 
xz 


Q 

o 


1 on 

lyu 


900 

ZUU 


QO 

OZ 


WSB 31 


77 
11 


J iDZDZ0.D-Z4Z4Zy 


1 R 


9R 91=1 RS 
ZD ZO.Do 


94 

Z4: 


94 
z^ 


9Q 
zy 


o 
o 


ZoU 


9S9 
Zc5Z 


4 Q 
4o 


LFAM 


4 




J iOzOzO.O-z4Z4oi 


1 fi 

J-U 


9fi 9fi "^S 


—94 

Z'i 


94 


■^1 

O J. 


n 
U 


QQfl 
OOU 


Q9Pi 
oZO 


OU 


LFAM 


5 


o 

U 


71(^9(^97 Q 9/19^1^0 
J IDZDZ / .o-Z4Zo0y 


1 R 
ID 


9R 97 7Q 
ZD z 1 . / y 


94 
— Z^ 


94 
Z^ 


n9 
uz 




z 


901 
ZUi 


90(^ 

zyo 


oO 








71(^9(^90 9/19Q17 
J iDZDZy.0-Z4Zoi / 


1 R 


9R 9Q 1^4 
ZD zy.o^ 


94 
— Z^ 


9"^ 
Zo 


1 7 
1 1 


Q 

o 


9/10 
Z4U 


97/1 
Z / 4 


QO 

oy 


LFAM 


7 




J ±0Z0oU.0-Z4ZZ00 




ZD OU.Ol 


94 
— Z'l 


99 
ZZ 


1^7 
( 


4 


iUDi 


1 OPi/l 
iU04 


QQ 
oy 


LFAM 


9 


llr 


J iOzOoi.O-Z4ioOa 




ZD Ol.Zo 


94 
— Z^ 


1 8 
lO 


'^9 
OZ 


Q 


9S20 
ZoU 


ooy 


Q8 
OO 


LFAM 


10 




71 fJ9fJQ9 7 9/199/1 P;* 

J iDZDoZ. <-Z4ZZ40 


1 R 
10 


9R '^9 71 
ZD OZ. / 1 


94 
— Z^ 


99 
ZZ 


4p; 

^0 


4 


900 

ZUU 


9PiO 

ZoU 


QO 

oy 








71 *^9*^QQ A 9/1191f^d 
J iDZDO<j.4-Z4iZiO 


16 


26 33.47 


-24 


12 


16 


u 


1 Q/in 
ie54U 


9090 
ZUZU 


QR 
oO 


SFAM 


12 




J iDZDoo.0-Z4Z44o 


16 


26 33.58 


-24 


24 


49 


n 
u 


91 O 

ziy 


907 
ZU / 


QO 

oy 








71 (^9(^Q/1 9 9/19Q9Q* 
J iDZDo4.Z-Z4ZoZo 


16 


26 34.20 


-24 


23 


28 


D 


71 onArA 

( iUUt^V ) 


1 1 /I Pi Pi 
11400 


QR 
oO 


LFAM 


11 




J162635.3-242405 


16 


26 35.34 


-24 


24 


05 


8 


758 


720 


36 


LFAM 


13 




J162639.0-243052 


16 


26 39.02 


-24 


30 


52 


6 


300 


288 


36 








J162642.3-242625 


16 


26 42.42 


-24 


26 


26 


8 


3364 


3560 


37 


LFAM 


15 


III 


J162643.8-241632* 


16 


26 43.77 


-24 


16 


33 


2 


2648 


2700 


50 


SFAM 


15 


III 


J162646.3-242001 


16 


26 46.32 


-24 


20 


02 





250 


270 


38 


LFAM 


17 




J162646.4-241935 


16 


26 46.40 


-24 


19 


35 


6 


223 


230 


38 








J162652.8-241953 


16 


26 52.89 


-24 


19 


53 


8 


288 


430 


44 








J162658. 7-242651 


16 


26 58.76 


-24 


26 


52 





1720 


2940 


51 


LFAM 


19 




J162700.0-242640* 


16 


26 59.99 


-24 


26 


40 


2 


55987 


57010 


51 


LFAM 


21 




J162700.0-242209 


16 


27 00.09 


-24 


22 


09 


4 


209 


340 


39 








J162702.1-241928 


16 


27 02.11 


-24 


19 


28 


6 


585 


549 


39 








J162702.3-242724 


16 


27 02.37 


-24 


27 


24 


4 


280 


360 


42 









''An asterisk following the source number indicates that additional notes follow below. VLA data (C- 
configuration) are from cleaned maps and fluxes are corrected for primary beam attenuation. Total flux 
is the average of the values measured using IMFIT (Gaussian fit) and TVSTAT (pixel summation) . A (v) 
following the total flux indicates that the flux is variable. RMS noise (±lcr) is measured in immediate 
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vicinity of source. Typical 6 cm radio position uncertainties are ±0.05 s in RA and ±0".7 in Dec. 

''FAM = Leous et al. 1991; SFAM = Stine et al. 1988; WLY = Wilking et al. (1989); WSB = Wilking 
et al. (1987) 

■^SED classification designations in order of preference are: Wilking et al. (2001, 1999); Luhman & Rieke 
(1999); Strom et al. (1995); Greene et al. (1994); Andre & Montmerle (1994). 

'^Outside the Chandra FOV. 
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Comments on individual sources: 

J162556. 1-243014- New radio detection. Probable counterpart is infrared source no. 11 in Table 2 of Wilking 
et al. (1989). 

J162603. 0-242335 - Emission at 6 cm is left circularly polarized (LCP) with a detection significance S/N = 16 
in Stokes V maps and a fractional circular polarization computed from peak Stokes V and I fluxes of 
= -0.059 (LCP). 

J162615. 5-243428 - Possible double radio source, but source structure is uncertain due to beam elongation 
at low source elevation. 

J162622.4-242252 - Peak flux increased by a factor of -5 during the time interval 0629 - 0931 UT. 
J162624. 0-241613 - New radio detection. Probable counterpart is Ha emission-line star no. 31 in Table 1 of 
Wilking et al. (1987), which is cross-Usted as Ehas-24 (EUas 1978). 

J162630. 5-242256 - The peak flux density at 3.6 cm is comparable to that at 6 cm, suggesting a flat spectral 

energy distribution and possible nonthermal radio emission. 
J162632.7-242245 - Marginal detection (S/N = 5.1). 

J162634. 1-242328 - Stokes I emission is extended. Emission at 6 cm is left circularly polarized (LCP) with a 
detection significance S/N = 9 in Stokes V maps and a fractional circular polarization tTc = —0.046 (LCP). 
Stokes V maps with Oph SI located near phase center give peak fluxes V(6 cm)= —330 (±38) /iJy and 
V(3.5 cm)= —390 (±30) /xJy, where the minus sign denotes LCP. The peak flux ratios V/I give fractional 
circular polarizations 7rc(6 cm)= —0.046 (LCP) and 7rc(3.5 cm)= —0.062 (LCP). The Stokes V map at 20 cm 
has an rms noise of 65 /uJy/beam, giving a 3a upper hmit V(20 cm) < 195 /iJy and 7rc(20 cm) < 0.0035. 
The 6-cm Stokes V flux given above is consistent with the value V(6 cm)= —280 (±50) /uJy obtained with 
the VLA on 1987 Aug 6 (Andre et al. 1988). 

J162643. 8-241632 - The infrared counterpart is oflfset by 3. "2 from the radio peak position making the 

identification uncertain. 

J162700. 0-242640 - Emission at 6 cm is right circularly polarized with a detection significance S/N = 34 
in Stokes V maps and a fractional circular polarization computed from peak Stokes V and I fluxes of 
TTc = +0.034. Corresponds to #20 (BZ6) in Table II of Stine et al. (1988), who argue that it is extragalactic. 



Table 4. X-ray sources in p Oph Core A with JHK Photometry 



GDS 


IR/Optical 


LFAM/ 


J 


H 


K 


log/x//K SED^ Av A^H 


logT 


log Lx 


Designation 


Designation 


SFAM 








(10^2 


(K) 


(ergs s"-*-) 



J162603.0- 242336 


DoAr 21 


SFAM 7 


8.01 


6.75 


6.16 


-1.75 


III 


6.6 


1 ic:+0.02 
— \).\JZ 


7. 


J162607.1-242724 


BKLT J162607-242725 




15.49 


12.50 


10.48 


-1.91 




19.5 


(3.47) 


7. 


J162607.6-242741 


BKLT J162607-242742 




14.73 


11.86 


10.27 


-1.97 




20.6 


(2.88) 


7. 


J162610.3-242054 


GSS 26 




14.97 


11.63 


9.38 


-2.71 


IIF 


22.6 


(4.17) 


7. 


J162615.8-241922 


SKS 1-7 




13.95 


11.44 


9.98 


-2.95 


III 


17.4 


(2.40) 


7. 


J162616.8-242223 


Oph S28 


LFAMp 1 


11.11 


9.23 


8.19 


-2.54 


II 


11.0 


(1.66) 


7. 


J162617.0-242021 


DoAr 24 




9.72 


8.67 


8.09 


-2.57 


II 


1.8 


0.24+»-!!^ 

— \.)A)I 


7. 


J162617. 1-241240 


BKLT J162617-241241 




13.11 


11.47 


10.65 


-2.86 






1.26+"-?; 

— (J. 7 1 


7. 


J162618.8-242819 


VSSG 1 




13.49 


10.76 


8.68 


-3.05 


II 


17.1 


2.40+"-!!^ 


7. 


J162621. 5-242600 


GY 5 




12.70 


11.57 


10.91 


-3.96 


III 


4.5 


(0.79) 


7. 


J162622.2- 242447 


SKS 3-12 




>17.00 


15.22 


12.97 


-3.05 






8.71+i\;?^ 

— ».(J5 


6. 


J162622.4-242252 


GSS 30#2 


LFAM 2 


14.22 


11.47 


9.60 


-1.64 


III 


lO.O'^ 


r, 00 + 0.02 
— U.Uz 


7. 


J162623.3-242059 


DoAr 24E A 




8.98 


7.53 


6.90 


-3.46 


II 


6.1 


(0.87) 


7. 


J162623.6-242439 


GY 21 


LFAM 3 


14.08 


11.63 


9.94 


-2.15 


IIF 


14.0 


2.63+0f2 


8. 


J162624. 0-242448 


Oph S2 




10.83 


8.68 


7.20 


-2.59 


II 


11.8 


3.47+;:-™ 

^ ' —0.01 


7. 


J162624. 0-241613 


Elias 24 




10.18 


8.23 


6.77 


-3.08 


II 


10.0 


(1.38) 


7. 


J162625. 2-242323 


GY 31 




>17.00 


> 15.50 


13.24 


-1.08 


III 


56.0 


(7.94) 


7. 


J162625. 2-242444 


GY 29 




15.92 


12.97 


10.86 


-3.53 


III 


19.1 


(3.47) 


7. 


J162625.7-241427 


BKLT J162625-241430 




13.04 


12.02 


11.48 


-2.74 






o.32+"';':f; 


7. 


J162627.8-242642 


GY 37 




14.25 


12.94 


11.99 


-3.43 


III 


6.3 


(1.05) 


7. 


J162628.5- 241540 


BKLT J162628-241543 




15.42 


12.43 


10.77 


-3.54 




21.7 


(2.88) 


7. 


J162629.7-241905 


SKS 1-19 


LFAM 8 


16.61 


13.24 


11.19 


-2.08 


III 


23.0 


(3.16) 


7. 


J162630.5-242256 


GY 51 


LFAM 9 


16.69 


13.46 


10.72 


-2.07 


IIF 


22.0 


(3.47) 


7. 


J162630.9-243106 


GY 54 




14.79 


12.36 


10.80 


-3.37 




16.6 


(2.19) 


7. 


J162631. 3-242530 


GY 59 




14.75 


12.89 


11.68 


-2.95 


III 


11.0 


(1.51) 


7. 


J162634.1-242328 


Oph SI B 




-9.00 


-9.00 


8.30 


-2.42 




10.0 


(1.82) 


7. 


J162636.8-241551 


BKLT ,1162636-241554 




12.13 


10.44 


9..38 


-3.12 




7.6 


1 OO+0.02 
-'-•■J°-0.02 


8. 


J162636.8-241900 


BKLT J162636-241902 




16.07 


13.93 


12.31 


-3.27 




14.0 


(2.63) 


7. 


J162637. 1-241601 


BKLT J162637-241602 




12.49 


11.37 


10.73 


-3.19 








7. 


J162640.4- 242714 


GY 91 




>17.00 


16.40 


12.51 


-1.71 






'■^^-3.44 


7. 


J162642.3-242625 


GY 101 


LFAM 15 


>17.00 


16.34 


12.58 


-2.33 


III 


50.0 


(6.61) 


7. 


J162642.5-242631 


GY 103 




>17.00 


15.82 


12.52 


-1.88 




48.5 


(6.61) 


7. 


J162642.8-242029 


GY 110 




10.71 


8.93 


8.00 


-2.83 


II 


8.5 


(1.51) 


7. 


J162642.9-242259 


GY 109 


LFAM 16 


15.37 


12.91 


11.39 


-3.56 




14.6 


(2.63) 


7. 



7.55 
7.45 
7.40 
7.55 
7.25 
7.40 
7.35 
7.15 
7.25 
7.50 
6.85 
7.50 
7.50 
8.35 
7.70 
7.60 
7.30 
7.00 
7.40 
7.15 
7.30 
7.35 
7.40 
7.40 
7.45 
7.40 
8.35 
7.60 
7.35 
7.45 
7.20 
7.35 



0.02 
0.02 
0.02 
-0.02 
0.02 
0.05 
0.30 
-0.05 
0.05 
0.10 
0.02 
-0.02 
0.02 
0.05 
0.30 
-0.15 
0.05 
0.02 
0.85 
-0.25 
1.50 
0.50 
0.05 
-0.02 
0.02 
0.05 
0.02 
-0.02 
0.02 
0.15 
0.02 
-0.05 
0.02 
0.10 
0.20 
-0.15 
0.50 
0.25 
0.40 
-0.35 
0.25 
0.25 
0.05 
-0.05 
0.10 
0.05 
+ 0.35 
-0.35 
0.20 
0.25 
0.02 
-0.02 
-0.02 
-0.02 
+ 0.75 
-0.30 
0.40 
O.fiO 
0.70 
-0.30 
-0.05 
-0.20 
-0.10 
-0.15 
+0.02 
-0.02 
+0.75 
-0.15 



31.82 
30.18 
30.20 



30.08Z 
29.22j 
29.87j 
27.67j 



29.79 



-0.02 
-0.02 
-0.02 
-0.02 
-0.07 
-0.02 

29.791;',;;? 

on ..c; + 0.02 

K).02 
-0.02 
-1 + 0.37 
^-0.37 
f0.02 
-O.Ofi 
7+0.12 
'-0.13 

29.73l|?;5 
30.721-- 

-0.03 
-0.02 

29.87i;;:!!^ 

on cv+O.ll 
30.67I(,„2 

29.i2i«:;?^ 

28.57l!5:i5 

28.05l5j;29 

28.52_(j 

29-79li',:o5 

29.96l!',-'i'*; 

28.53l«;- 

28.49ii;:;5:^ 

f0.02 
*-0.02 
f0.02 
-0.02 

D+i:o7 



30.44Z 
29.06j 



28.38_m^ 
29 75+''-'^'^ 
29 45+0-™ 

29.67l°f5' 
30.061H2 

28.14l«:i« 



Table 4 — Continued 



GDS 


IR/Optical 


LFAM/ 


J 


H 


K 


log /x//k 


SED* 


Av 




logT 


logLx 


Designation 


Designation 


SEAM 














(1022cm-2) 


(K) 


(ergs s""') 


J162643.8-241632 


VSSG 11 


SEAM 15 


13.12 


10.85 


9.58 


-2.57 


III 


14.0 


(2.00) 


7 nc + O.OS 


on 04+0-09 


J162645.0- 242307 


GY 116 


LFAMp 3 


13.14 


10.50 


8.88 


-3.58 


II 


17.0 


(3.16) 




28-98±o.l2 


J162646.8-241907 


BKLT J162646-241910 




14.05 


11.67 


10.34 


-3.21 






^•°^-0.91 


7 Qc+o.ss 
' •■-''-'-0.25 


OO oo+0!49 


J162648.3-242834 


BKLT J162648-242836 




> 17.00 


15.73 


12.48 


-2.40 








8.2oi;:jg 


28.56lO:J« 


J162648.6-242839 


GY 128 




> 17.00 


14.05 


10.99 


-3.12 




38.6 


(7.24) 


7 qn+'^-^S 
'•"^'J-0.10 


29.13t°io 


J162649.2-242002 


GY 135 


LFAM 18 


12.15 


9.82 


8.62 


-2.56 


III 


16.0 


^•-■-^-0.02 


7 QK+0.02 


on 2^ + 04 


J162651. 9-243039 


GY 144 




>17.00 


15.75 


13.46 


-2.42 




26.8 


(3.80) 


7 cc+O.SO 
' .OO_0.40 

7 2n+°'^-'' 

' •^'-'-().20 


28.43l°;5| 


J162653.4-243236 


GY 146 




>17.00 


>15.50 


13.18 


-2.15 




43.5 


(6.03) 


29.35l°;i 


J162654.4-242620 


GY 153 




14.72 


11.70 


9.88 


-2.53 


III 


19.7 


(2.63) 




29.87lO:J2 


J162655.0-242229 


GY 156 




15.04 


11.97 


10.19 


-2.55 


III 


22.5 


(3.47) 


7 oQ+o.io 


29.79+„ 


J162658.3-242130 


GY 171 




16.12 


13.38 


11.74 


-2.65 




19.9 


(3.47) 


7 4c;+0-45 
'•*3-0.10 


28.9410:09 


J lDzOOo.b-Z4ioo4 


rSrvlji J lDZDOo-z41ooD 




1 C /I 

15.47 


13.06 


11.54 


-3.26 




14.2 


(2.63) 


7 r,c,+O.Xy 

"•20_„.io 


ZO.oU Q 3y 


J162700.5-241624 


BKLT J162700-241627 




13.46 


11.86 


10.95 


-3.15 






U.OU_Q 43 




28.44lo;« 


J162703.0-242614 


GY 188 




> 17.00 


>15.50 


12.52 


-2.21 




20.5 


(3.80) 


8 H5+° °2 


no 00 + 0.03 
2o-o3_Q 02 


J162704. 1-242828 


GY 192 




16.22 


13.07 


10.81 


-3.01 




20.9 


(3.80) 


7 4C-I-0.55 
'•^o_0.10 


28.971°°! 


J162704.6-242715 


GY 195 




16.67 


13.55 


11.30 


-0.75 




20.6 


3.80l°;°l 


7 60+°-°^ 


30.941°;°^ 


J162705.2- 242007 


GY 199 




12.48 


10.41 


9.27 


-3.05 






2.ool°;« 


7 40+0-35 


29 44+0-13 
^"•^^-0.22 


J162706.0-242618 


GY 203 




>17.00 


13.57 


11.05 


-2.18 




36.5 


(5.01) 


7 oQ+0.05 


29.931°:?^ 


J162710.3-241912 


Elias 30 A 




8.56 


7.33 


6.30 


-3.93 


IID 


3.5 


U.OO_Q gg 


7 20+"-"-''' 

'■^'-'-().02 


29.62lo-0| 


J162710.3-241918 


Elias 30 B 




11.27 


10.00 


9.25 


-3.64 




3.5 


(0.50) 


7 3Q+o.i() 


28.661°;°! 


J162713.7-241816 


VSSG 24 




12.30 


10.27 


9.32 


-3.81 




10.7 


(1.82) 


7.30l°:?g 


28.721°:!^ 


J162722.9-241800 


VSSG 22 




13.27 


10.79 


9.41 


-2.52 


III 


16.0 


(2.88) 


7 on-H0.05 
' "^"-0.02 


30.09lo 07 



Note. — The quoted errors arc internal errors only based on the non-parametric method described in Appendix B. Additional errors due to systematic 
effects such as distance and Ay uncertainties are not included. 

Note. — All JHK photometry from Barsony ct al. (1997) except Oph 81 B, K=8.3 (Simon et al. 1995). Lx is the unabsorbed 0.5-7.0 keV luminosity 
assuming d = 165 pc. Parentheses indicate Nn estimated from Ay assuming 1.57 X lO'^icm"^ per mag (Vuong et al. 2003). 

■^D = double peaked spectrum; F = flat spectrum. 

''Brown dwarf candidate. 



°Near-IR measurements give Ay f» 27 mag (Casteiaz et al. 1985) 
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Table 5. X-ray Upper Limits for Class O-III Objects in p Oph 



Name 


Class 


log Nif 
(cm-2) 


log Lx (0.5 - 7 keV) 
(ergs s"^) 


LFAM 5 





22. IS-'' 


<27.81 


GSS 30-1 


I 


22.78'' 


<28.32 


GSS 30-3 


I 


22.48=' 


<28.05 


GY 81 


Il(flat) 


22.46 


<28.00 


GY 84 


IIP 


22.40 


<27.71 


GY 154 


II 


22.59 


<28.27 


SKS 1-11 


IIP 


22.48 


<27.92 


''Column density is an assumed value that is typical of 
sources in p Oph core-A. 

''Based on Ay = 38 (Castelaz et al. 1985) and the con- 



version of Vuong ct al. (2003). 
'^Classification uncertain. 
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Table 6. Spectral Fit Results for Oph Sl^ 



Parameter 


PL 




IT 




2T 




Nh (10^2 cm-2) 


2.16 [1.99 - 


2.36] 


1.55 [1.48 - 


1.68] 


1.90 [1.50 - 


2.45] 


kTi (keV) 






2.38 [2.17 - 


2.60] 


0.46 [0.29 - 


1.95] 


kT2 (kcV) 










2.32 [2.03 - 


2.69] 


Photon index 


3.00 [2.84 - 


3.19] 










Fe 






0.10 [0.00 - 


0.33] 


0.24 [0.00 - 


0.57] 


EM2/EM1 










1.79 




xVdof 


147.0/135 




141.2/134 




137.3/132 






1.09 




1.05 




1.04 




ix (IQ--^^ ergs cm-^ s-^) 


3.00 




2.99 




3.01 




log Lx (ergs s'^) 


30.75 




30.37 




30.59 





''PL = power-law model. IT and 2T denote VAPEC optically thin plasma models 
with one and two temperature components. Fe abundance is relative to solar and abun- 
dances for all other elements were held fixed at their solar values (Anders & Grevesse 
1989). ix is the absorbed flux (0.5 - 7 keV) measured from model fit. Unabsorbed 
luminosity Lx (0.5 - 7 keV) assumes d = 165 pc. Brackets enclose 90% confidence 
ranges. Spectra are corrected for charge transfer ineflaciency (CTI) using CIAO v. 2.3 
and CALDB v. 2.20 and rebinned to a minimum of 15 counts per bin. 
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Table 7. MARX VAPEC Simulations: Nu unknown 



Input Model Output Parameters 



lop- AAxj 
lug i V H 


loa-T 
lug -t 


FCF 


N 


lug JVJJ 


; 


logT 




FCF 


FCF 68% c.l. 


(cm j 






[CZS ) 


(cm 


(K) 








21.50 


7.0 


2.04 


10 


21.26 


0.65 


6.92 


0.33 


1.68 


1.06 


53.28 


21.50 


7.4 


1.47 


10 


21.60 


0.58 


7.32 


0.51 


1.66 


1.12 


14.74 


21.50 


7.8 


1.28 


10 


21.62 


0.60 


7.65 


0.56 


1.37 


1.10 


6.36 


22.18 


7.0 


7.62 


10 


22.19 


0.41 


6.92 


0.35 


11.26 


2.42 


286.65 


22.18 


7.4 


2.53 


10 


22.27 


0.36 


7.36 


0.46 


3.64 


1.52 


24.16 


22.18 


7.8 


1.79 


10 


22.31 


0.35 


7.70 


0.51 


2.04 


1.51 


10.08 


22.90 


7.0 


51.41 


10 


22.92 


0.26 


6.99 


0.35 


87.53 


6.52 


4493.47 


22.90 


7.4 


6.33 


10 


22.97 


0.25 


7.45 


0.46 


7.22 


2.57 


207.73 


22.90 


7.8 


3.31 


10 


23.05 


0.23 


7.73 


0.44 


4.52 


2.57 


26.98 


21.50 


7.0 


2.04 


25 


21.29 


0.50 


6.97 


0.21 


1.69 


1.13 


4.11 


21.50 


7.4 


1.47 


25 


21.45 


0.51 


7.37 


0.31 


1.43 


1.12 


4.52 


21.50 


7.8 


1.28 


25 


21.56 


0.39 


7.79 


0.41 


1.33 


1.14 


1.96 


22.18 


7.0 


7.62 


25 


22.15 


0.21 


7.00 


0.22 


7.62 


3.36 


19.39 


22.18 


7.4 


2.53 


25 


22.17 


0.22 


7.46 


0.33 


2.47 


1.60 


5.89 


22.18 


7.8 


1.79 


25 


22.21 


0.18 


7.81 


0.37 


1.94 


1.46 


2.86 


22.90 


7.0 


51.41 


25 


22.91 


0.17 


6.98 


0.20 


58.52 


14.20 


541.86 


22.90 


7.4 


6.33 


25 


22.90 


0.22 


7.49 


0.39 


5.27 


2.51 


58.16 


22.90 


7.8 


3.31 


25 


23.01 


0.19 


7.73 


0.40 


4.03 


2.57 


18.50 


21.50 


7.0 


2.04 


60 


21.33 


0.42 


6.98 


0.14 


1.86 


1.18 


3.71 


21.50 


7.4 


1.47 


60 


21.49 


0.30 


7.40 


0.22 


1.43 


1.22 


1.92 


21.50 


7.8 


1.28 


60 


21.48 


0.30 


7.81 


0.28 


1.30 


1.15 


1.53 


22.18 


7.0 


7.62 


60 


22.16 


0.11 


7.00 


0.12 


7.62 


4.20 


13.98 


22.18 


7.4 


2.53 


60 


22.16 


0.13 


7.46 


0.19 


2.32 


1.87 


3.64 


22.18 


7.8 


1.79 


60 


22.19 


0.13 


7.86 


0.31 


1.81 


1.52 


2.40 


22.90 


7.0 


51.41 


60 


22.88 


0.12 


7.01 


0.14 


45.27 


13.04 


199.65 


22.90 


7.4 


6.33 


60 


22.87 


0.16 


7.49 


0.29 


5.39 


2.81 


16.81 


22.90 


7.8 


3.31 


60 


22.96 


0.13 


7.82 


0.32 


3.63 


2.71 


6.33 


21.50 


7.0 


2.04 


160 


21.40 


0.27 


7.00 


0.08 


1.85 


1.42 


2.60 


21.50 


7.4 


1.47 


160 


21.48 


0.16 


7.42 


0.11 


1.43 


1.30 


1.63 


21.50 


7.8 


1.28 


160 


21.48 


0.15 


7.86 


0.19 


1.28 


1.17 


1.38 


22.18 


7.0 


7.62 


160 


22.17 


0.06 


7.00 


0.08 


7.62 


4.86 


12.56 


22.18 


7.4 


2.53 


160 


22.16 


0.08 


7.43 


0.13 


2.43 


2.01 


3.12 


22.18 


7.8 


1.79 


160 


22.16 


0.08 


7.90 


0.23 


1.71 


1.55 


2.04 


22.90 


7.0 


51.41 


160 


22.88 


0.08 


7.01 


0.10 


45.27 


21.09 


142.62 


22.90 


7.4 


6.33 


160 


22.87 


0.10 


7.46 


0.15 


5.27 


3.77 


8.79 


22.90 


7.8 


3.31 


160 


22.93 


0.07 


7.85 


0.23 


3.31 


2.71 


4.39 


21.50 


7.0 


2.04 


400 


21.44 


0.15 


7.01 


0.04 


1.94 


1.69 


2.29 


21.50 


7.4 


1.47 


400 


21.45 


0.10 


7.43 


0.07 


1.40 


1.30 


1.55 


21.50 


7.8 


1.28 


400 


21.47 


0.11 


7.85 


0.16 


1.28 


1.19 


1.33 


22.18 


7.0 


7.62 


400 


22.17 


0.03 


7.01 


0.05 


7.62 


5.85 


9.25 


22.18 


7.4 


2.53 


400 


22.16 


0.06 


7.42 


0.08 


2.43 


2.09 


2.87 


22.18 


7.8 


1.79 


400 


22.15 


0.06 


7.90 


0.18 


1.71 


1.55 


1.89 


22.90 


7.0 


51.41 


400 


22.88 


0.06 


7.02 


0.08 


51.41 


21.09 


70.98 
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Table 7 — Continued 



Input Model Output Parameters 



logATn logT FCF N logATn cj^ logT cTjt FCF FCF 68% c.l. 
(cm-2) (K) (cts) (cm-2) (K) 



22.90 7.4 6.33 400 22.87 0.08 7.46 0.11 4.98 3.90 8.54 

22.90 7.8 3.31 400 22.92 0.06 7.86 0.21 3.31 2.73 4.03 
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Table 8. MARX VAPEC Simulations: 10% A^h error 



Input Model Output Parameters 



log iViJ 


log 1 


FOF 




log YVn 


" ) 


logT 




FCF 


FCF 68% c.l. 


((-m - ) 


(iV ) 




Is) 


(t'lll 










21.50 


7.0 


2.04 


10 


21.49 


0.04 


7.01 


0.17 


2.04 


1.68 


2.34 


21.50 


7.4 


1.47 


10 


21.50 


0.04 


7.48 


0.31 


1.47 


1.30 


1.66 


21.50 


7.8 


1.28 


10 


21.51 


0.04 


7.80 


0.39 


1.31 


1.22 


1.47 


22.18 


7.0 


7.62 


10 


22.19 


0.04 


6.98 


0.22 


8.35 


4.00 


18.88 


22.18 


7.4 


2.53 


10 


22.19 


0.04 


7.47 


0.30 


2.63 


1.89 


4.00 


22.18 


7.8 


1.79 


10 


22.19 


0.03 


7.81 


0.36 


1.80 


1.62 


2.47 


22.90 


7.0 


51.41 


10 


22.90 


0.04 


7.00 


0.15 


58.52 


17.19 


169.44 


22.90 


7.4 


6.33 


10 


22.90 


0.04 


7.54 


0.37 


5.94 


2.95 


12.83 


22.90 


7.8 


3.31 


10 


22.92 


0.03 


7.88 


0.39 


3.31 


2.79 


5.94 


21.50 


7.0 


2.04 


25 


21.49 


0.03 


7.01 


0.08 


2.04 


1.85 


2.19 


21.50 


7.4 


1.47 


25 


21.49 


0.04 


7.45 


0.21 


1.43 


1.33 


1.59 


21.50 


7.8 


1.28 


25 


21.51 


0.04 


7.84 


0.30 


1.29 


1.22 


1.38 


22.18 


7.0 


7.62 


25 


22.18 


0.04 


6.96 


0.14 


8.35 


4.86 


15.21 


22.18 


7.4 


2.53 


25 


22.18 


0.04 


7.43 


0.17 


2.47 


2.07 


3.36 


22.18 


7.8 


1.79 


25 


22.18 


0.04 


7.90 


0.29 


1.75 


1.61 


2.07 


22.90 


7.0 


51.41 


25 


22.90 


0.04 


6.99 


0.09 


58.52 


29.05 


122.10 


22.90 


7.4 


6.33 


25 


22.90 


0.04 


7.43 


0.18 


6.75 


3.96 


10.16 


22.90 


7.8 


3.31 


25 


22.91 


0.03 


7.92 


0.33 


3.22 


2.79 


4.39 


21.50 


7.0 


2.04 


60 


21.50 


0.03 


7.00 


0.05 


2.04 


1.94 


2.19 


21.50 


7.4 


1.47 


60 


21.49 


0.04 


7.41 


0.12 


1.47 


1.37 


1.56 


21.50 


7.8 


1.28 


60 


21.50 


0.04 


7.83 


0.24 


1.30 


1.22 


1.36 


22.18 


7.0 


7.62 


60 


22.18 


0.04 


7.00 


0.08 


7.62 


5.65 


10.20 


22.18 


7.4 


2.53 


60 


22.18 


0.04 


7.40 


0.10 


2.63 


2.21 


3.12 


22.18 


7.8 


1.79 


60 


22.18 


0.04 


7.87 


0.23 


1.75 


1 R9 




22.90 


7.0 


51.41 


60 


22.89 


0.03 


7.00 


0.08 


51.41 


32.06 


90.28 


22.90 


7.4 


6.33 


60 


22.90 


0.03 


7.43 


0.12 


5.94 


4.31 


8.16 


22.90 


7.8 


3.31 


60 


22.90 


0.03 


7.96 


0.28 


3.16 


2.73 


3.65 


21.50 


7.0 


2.04 


160 


21.49 


0.03 


7.00 


0.03 


2.04 


1.94 


2.19 


21.50 


7.4 


1.47 


160 


21.49 


0.03 


7.41 


0.06 


1.43 


1.40 


1.50 


21.50 


7.8 


1.28 


160 


21.49 


0.03 


7.82 


0.14 


1.28 


1.24 


1.32 


22.18 


7.0 


7.62 


160 


22.18 


0.03 


6.99 


0.06 


7.62 


5.85 


10.20 


22.18 


7.4 


2.53 


160 


22.18 


0.03 


7.40 


0.07 


2.53 


2.29 


2.87 


22.18 


7.8 


1.79 


160 


22.17 


0.03 


7.88 


0.18 


1.73 


1.62 


1.95 


22.90 


7.0 


51.41 


160 


22.89 


0.03 


7.01 


0.05 


45.27 


32.06 


70.98 


22.90 


7.4 


6.33 


160 


22.89 


0.03 


7.43 


0.08 


5.94 


4.98 


7.41 


22.90 


7.8 


3.31 


160 


22.89 


0.03 


7.98 


0.24 


3.16 


2.71 


3.31 


21.50 


7.0 


2.04 


400 


21.50 


0.03 


6.99 


0.02 


2.04 


1.94 


2.16 


21.50 


7.4 


1.47 


400 


21.48 


0.03 


7.40 


0.04 


1.43 


1.40 


1.51 


21.50 


7.8 


1.28 


400 


21.49 


0.03 


7.82 


0.10 


1.28 


1.24 


1.30 


22.18 


7.0 


7.62 


400 


22.17 


0.03 


7.00 


0.05 


7.62 


5.85 


9.25 


22.18 


7.4 


2.53 


400 


22.17 


0.03 


7.41 


0.05 


2.43 


2.29 


2.87 


22.18 


7.8 


1.79 


400 


22.16 


0.03 


7.88 


0.10 


1.68 


1.63 


1.84 


22.90 


7.0 


51.41 


400 


22.89 


0.03 


7.01 


0.03 


51.41 


32.06 


51.41 



-50- 



Table 8 — Continued 



Input Model Output Parameters 



logATn logT FCF N logATn cj^ logT cTjt FCF FCF 68% c.l. 
(cm-2) (K) (cts) (cm-2) (K) 



22.90 7.4 6.33 400 22.89 0.02 7.42 0.06 6.33 5.27 6.93 

22.90 7.8 3.31 400 22.89 0.03 7.98 0.23 3.16 2.62 3.31 
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Fig. 1. — Pull-resolution grayscale ACIS-I image of p Ophiuchus cloud A (L1688) with a 
logarithmic stretch. Pixel size is 0.492". Sources listed in Table 1 are indicated with open 
circles and ellipses. The ACIS-I aimpoint is the magnetic B4 star Oph SI at 16'^26™34® — 
24°23^5. The weak-hned T Tauri star DoAr 21 at far right (16^26°'03" - 24°23;6) is in the 
decay stage of a large flare throughout the ~ 96 ks observation. The trailed image near 
DoAr 21 is caused by the 40/isec row-to-row charge transfer time. 

Pig. 2. — Light curves of /x, logT and Lx (assuming d = 165 pc) for the optically bright 
WTTS DoAr 21=J162603.0-242336 (KO, class III). The bin size is varied so as to include 
approximately the same number of counts per bin (and similarly for Pigurcs 3-8). The 
flare decay is punctuated by smaller heating events with coronal temperatures varying from 
^ 28 - 50 MK. 

Pig. 3.— Light curves of /x, logT and Lx for the CTTS Oph S2=J162624.0-242448 (class 
II). Oph S2 underwent a very rapid temperature spike (T > 100 MK) at UT 19h with the 
luminosity peak occuring nearly two hours later. The luminosity decayed with an e-folding 
time of approximately 6.5h and a steady decrease in temperature from > 100 MK to 25 MK. 

Pig. 4. — Light curves of /x, logT and Lx for the B4-|-K-type binary Oph S1=J162634.1- 
242328. The presence of hot coronal plasma (T 19 MK) and flare-like variability is typical 
of Class II or Class III TTSs. It is likely that most of Oph Si's X-ray emission is produced 
by the low- mass K-type companion discovered via lunar occultation by Simon et al. (1995), 
not the magnetic B star of Andre et al. (1988). 

Pig. 5.— Light curves of /x, logT and Lx for GY 195=J162704.6-242715. This imclassiflcd 
star underwent the largest flare-amplitude variation seen in our observation: a twenty-fold 
increase in X-ray luminosity over Ih followed by a 2h decay. The temperature spikes (which 
exceed 50 MK) after 13h UT do not suggest rotational modulation of the inital large flare; 
rather they suggest two to three subsequent reheating events. 

Pig. 6.— Light curves of /x, logT and Lx for the CTTS Oph S28=J162616.8-242223 (K6, 
class II). The initial flare has rise and decay times of approximately 2h and 6h, respectively. 

Pig. 7.— Light curves of /x, logT and Lx for the CTTS DoAr 24=J162617.0-242021 (class 
II). This CTTS hght curve shows moderate- intensity, relatively short-duration flares through- 
out. 



Pig. 8.— Light curves of /x, logT and Lx for the WTTS GSS 30 IRS-2 = J162622.4-242252 
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(K8, class III). Typical WTTS light curves like this one are similar to those of the CTTS 
DoAr 24 and the unclassified K-type star Oph SI B. The persistent increases in mean energy 
may be the result of nearly continuous reheating and the relatively short rise and decay times 
suggest moderate-sized magnetic loops. 

Fig. 9. — (a) Kolmogorov-Smirnov (KS) variabhhty statistic versus mean energy ((E)) for 
flat or Class II sources (squares) , Class III sources (triangles) , stars with no SED classification 
(circles), and Chandra sources without NIR counterparts (crosses). Brown-dwarf candidates 
(BDC) are indicated with gray symbols. Chandra sources without NIR counterparts are 
typically very hard ((E) > 3.5. marked by vertical dotted line) and not variable (KS* < 1.0) 
strongly suggesting that these arc highly absorbed cxtragalactic sources, (b) count rate 
versus (E), and (c) count rate versus KS statistic for all Chandra sources in Table 1. A 
dashed fine at KS = 1 in (a) and (c) marks the statistical variabihty threshold. 

Fig. 10. — Near-infrared color-magnitude diagram for stars in the ACIS-1 field of view. 
Only stars with measured J, H. and K arc shown. A right arrow denotes objects whose J 
magnitudes are lower limit. Stars detected by Chandra are indicated with filled symbols 
and undetected stars with open symbols. Class I (diamonds), fiat or Class II (squares), and 
Class III (triangles) designations are based on the spectral energy distributions of Luhman 
& Rieke (1999), Wilking et al. (2001), Andre & Montmerle (1994), and Strom et al. (1995). 
The known Class source LFAM 5 was not detected at K-band or with Chandra, and the 
undetected class I source GSS 30-3 (not shown) has no measured J or H magnitudes. The 
dotted polygon encloses the region occupied by several dozen undetected near-IR sources in 
the Chandra FoV which lack published SED classifications. 

Fig. 11. — Near- infrared color-color diagram for stars in the ACIS-I FOV. Only stars with 
measured J, H, and K are shown. Symbol definitions are the same as in the previous figure. 
Right and up arrows denote objects whose J and H magnitudes are lower limits, respectively. 
The Ay — 10 reddening vector shown is typical for sources in p Oph core A. The dotted 
polygon encloses the region occupied by several dozen undetected near-IR sources in the 
Chandra FoV which lack published SED classifications. The unreddened ZAMS is shown as 
a dash-dot line at lower left. The two sloping dashed lines mark the approximate reddening 
band for normally reddened MO V (left) and AO V (right) stars using data from Bessell & 
Brett (1988) and Rieke & Lebofsky (1985). 

Fig. 12. — X-ray to IR flux ratio for all sources with measured K-band magnitude. The 
dotted line represents the total data while the solid line indicates bona flde YSOs with SED 
classiflcations. 

Fig. 13. — Kaplan-Meier estimator for 14 CTTS (12 detections and 2 upper limits) and 17 
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WTTS (15 detections and 2 upper limits) in the p Oph core-A ACIS-I field. The generahzed 
Wilcoxon test gives a probability p — 0.53 that the luminosities are drawn from the same 
parent distribution. 

Fig. 14. — Lq cm radio luminosity versus Lx (0.5 - 7 keV) for class II (CTTS) and class III 
(WTTS) sources detected simultaneously with Chandra and the VLA. We assume a distance 
of 165 pc. The unusual binary Oph SI and three class O/I sources are also labeled. Only 
X-ray upper limits were obtained for the class O/I sources. GSS30-IRS1 was not detected 
by Chandra or the VLA. The quiescent X-ray and radio luminosities of coronal active stars 
(main-sequence and evolved giants) generally obey the relation log Lx < log Lj. (Giidcl, M. 
& Benz 1993) and are located to the left of the dashed line. The PMS stars shown here 
also satisfy this inequality. The five X-ray brightest VLA sources underwent flares during 
the Chandra observation. Their Lx ranges are indicated with horizontal lines. The Class III 
sources appear to be brighter at 6 cm than other types of YSOs. The location of Oph SI in 
this diagram suggests that both the X-ray and radio emission are produced on the low-mass 
companion and that Oph SI B is a Class III source. 

Fig. 15. — Background-subtracted ACIS-I spectrum of the binary system Oph SI (B4 + K). 
The spectrum is rebinned to a minimum of 15 counts per bin. Solid line shows a best-fit 2T 
VAPEC optically thin plasma model with a variable Fe abundance. Best-fit values are given 
in Table 6. 

Fig. 16. — Background-subtracted ACIS-I light curve of the brown dwarf candidate GY 31. 
Nearly all the flare emission is in the 2-7 keV band. Although the flrst flare peak was not 
observed, fiare decay times appear to be less than five hours, typical of moderate fiares on 
other TTS. 

Fig. 17. — Rebinned background-subtracted ACIS-I spectrum of the brown dwarf candidate 
GY 31. Sohd line is a best-fit IT VAPEC model with log A'h = 22.76 cm'^, kT = 2.3 keV, 
Fe abundance Fc = 0.3 solar, and x^/dof = 19/26. 

Fig. 18. — logT versus logA^n contour plots for two fake sets of photon events: a source 
with 25 counts (upper panel), another with 60 counts (lower panel). The input values were 
logA^H = 22.18 and logT = 7.40 (indicated with crosses). The best-fit values from the 
non-parametric method described in Appendix B are indicated with dots. Also shown are 
the 68%, 90% and 95% confidence contours for the Cramer- von Mises statistic. Multiple 
fake data sets were used to estimate the standard deviations of the best-fit logA^n and 
logT = 7.40, indicated as a dashed box. This analysis shows that the non-parametric 
method can be used to estimate column density, temperature, and Lx with as few as 25 or 
60 counts. See Appendix B and Tables 7 and 8 for details. 
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